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A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly 
collisional plasma (e.g., drift-wave turbulence driven by ion temperature gradients) is 
proposed, with account taken both of the nonlinear advection of the perturbed particle 
distribution by fluctuating E x B flows and of its phase mixing, which is caused by 
the streaming of the particles along the mean magnetic field and, in a linear problem, 
would lead to Landau damping. It is found that it is possible to construct a consistent 
theory in which very little free energy leaks into high velocity moments of the distribution 
function, rendering the turbulent cascade in the energetically relevant part of the wave- 
number space essentially fluid-like. The velocity-space spectra of free energy expressed 
in terms of Hermite-moment orders are steep power laws and so the free-energy content 
of the phase space does not diverge at infinitesimal collisionality (while it does for a 
linear problem); collisional heating due to long-wavelength perturbations vanishes in this 
limit (also in contrast with the linear problem, in which it occurs at the finite rate 
equal to the Landau-damping rate). The ability of the free energy to stay in the low 
velocity moments of the distribution function is facilitated by the “anti-phase-mixing” 
effect, whose presence in the nonlinear system is due to the stochastic version of the 
plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing 
perturbations). The partitioning of the wave-number space between the (energetically 
dominant) region where this is the case and the region where linear phase mixing wins its 
competition with nonlinear advection is governed by the “critical balance” between linear 
and nonlinear timescales (which for high Hermite moments splits into two thresholds, one 
demarcating the wave-number region where phase mixing predominates, the other where 
plasma echo does). 


1. Introduction 

Turbulence is a process whereby energy injected into a system (via some mechanism 
usually associated with the system being out of equilibrium) is transferred nonlinearly— 
and therefore leading to chaotic and multiscale states—from the scale(s) at which it is 


f Email: alex.schekochihin@physics.ox.ac.uk 


2 


A. A. Schekochihin et al. 


injected to much smaller scales at which it is thermalised through microphysical dissipa¬ 
tion channels available in the system. The system is forced to seek ways of transferring 
energy across a range of scales because the injection and dissipation physics are usually 
unrelated to each other and operate at disparate scales. It is the bridging of the gap be¬ 
tween these scales that brings about turbulent cascades, broad-range power-law spectra, 
and so on. In fluid systems, however varied and multi-physics they a re, most turbulence 
the ories are basically ex tensions and generalisations of the ideas of iRichardson (1922) 
and iKolmogorovI (1194161 ) of a local-in-scale cascade maintaining a constant flux of energ:v 
away from the injection and towards the dissipation scales fe.g.. IZakharov et al. 


e nergy 

]ll992l: 


Davidsonll2013f) . This type of thinking has been tremendously successful in making sense 


of experimental and numerical evidence in both fluids and plasmas. 

In plasmas, however, a straightforward application of such “fluid” thinking to any 
physical regime that is not collisionally dominated skirts over the obvious complication 
that the kinetic phase space includes the particle velocities as well as their positions, 
and the (free) energy is generally free to travel across this entire 6D space. Its ability— 
and propensity—to do so is, in fact, manifest in what is probably the most important 


m portant 
3 (Il946h 


phenomenon that makes plasmas conceptually different from fluids—the iLandar 
damping of electromagnetic perturbations in a collisionless plasma. Viewed in energy 
terms, it involves the transfer of free energy from electromagnetic perturbations into 
perturbations of the particle distribution function, which develops ever finer structure in 
velocity space (“phase mixing”) until this transfer (which looks like damping if one only 
tracks the electromagnetic fields) is made irreversible by coarse graining of the velocity- 
space structure. The physical agent of this coarse graini ng is coll i sions, even if they would 
appear to be infinitesimally small. Mathematically, the iLandaiJ ( 19361) collision operator 
is a diffusion operator in velocity space and so even small collision frequencies are enough 
to thermalise any amount of energy, provided sufficiently large velocity-space gradients 
develop. 

In a linear plasma system, Landau damping, or, more generally, phase mixing, is the 
only available thermalisation route. It provides an adequate mechanism to process any 
injected free energy at any fixed wave number (since the process is linear, energy will 
stay in the wave number into which it is injected; there is no coupling), leading to a 
finite effective damping rate and filling up the phase space with free energy. If one uses a 
Hermite decomposition to quantify “scales” in velocity space, one finds that, in a steady- 
state system continuously pumped via low Hermite moments and dissipating free energy 
via high ones, the free energy will accumulate in phase space to a level that diverges if 
the collisionality is taken to zer o; the collisional hea ting rate in this limit is finite and 
equal to the phase-mixing rate ( Kanekar et al. [iol^. How does this mechanism coexist 
and compete with the refinement of spatial scales caused by coupling between scales—a 
well-nigh inevitable consequence of nonlinearity? 

In this paper, we address this question using a simple archetypal example of plasma 
turbulence—electrostatic turbulence in a drift-kinetic plasma. We will describe this ex¬ 
ample in section [21 along with all the relevant preliminaries: the concept of free energy, 
the Hermite deco mposition, and the e xisting Kolmogorov-style “fluid” turbulence theory 
for this problem (jBarnes et aLl 120111 ). In section |3l we will introduce the phase-space 
formalism that explicitly separates the phase-mixing and the “anti-phase-mixing” per¬ 
turbations (the latter activated by the plasma echo effect), both of which turn out to be 
inevitable in a nonlinear system, and provides a useful starting point for a substantive 
theoretical treatment of phase-space turbulence. In section |4l a phenomenological scaling 
theory of this turbulence will be proposed. While we will describe in detail how free 
energy and its fluxes are distributed in the inertial range—leading to some interesting 
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and testable scalings—the main conclusion will be that phase mixing is quite heavily 
suppressed in a turbulent system. Section [5] is devoted to summarising this and other 
findings and to discussing their implications, as well as future directions of travel. A 
reader only interested in a digest can skip to this section now. 


2. Preliminaries 

This section contains a rather extended tutorial on a number of topics constituting 
elementary but necessary background to what will follow. Readers who are sufficiently 
steeped in these matters can skim through this section and then dedicate themselves 
more seriously to sections [3] and S] (where references to relevant parts of section [3] will be 
supplied). 


2.1. Prototypical kinetic problem 

We consider a plasma near Maxwellian equilibrium, in which case the distribution func¬ 
tion for particles of species s can be expressed as 

fs=FMs+5fs, ( 2 . 1 ) 

where Tms is a Maxwellian distribution and Sfg a small perturbation. 

We assume this plasma to be in a uniform strong magnetic field B = Bz (z is the unit 
vector in the direction of this field, designated the z axis). We consider low-frequency 
perturbations, which will be highly anisotropic with respect to the field: 

uj < Ps, fc|| < fc_L, (2.2) 

where fls is the Larmor frequency. 

We assume these perturbations to be electrostatic, viz., 

(5E = - V(/), (5B = 0, (2.3) 


where (f) is the scalar potential. We use Gaussian electromagnetic units. 
We consider only long wavelengths, 

k±Ps < 1 , 


(2.4) 


where ps is the Larmor radius. 

Finally, we assume a Boltzmann electron response (which arises via expansion of the 
electron drift-kinetic equation in the electron-to-ion mass ratio) U 


e(l) 6ne Sui ^ f ,3 c n 

7^ = -=- = — i vSfi, 

If. Ue Ui Ui J 


(2.5) 


where e is the electron charge, Tg and are the equilibrium temperatures and number 
densities, respectively, and Shs are density perturbations (s = e for electrons, s = i for 
ions). The second equality in equation (12.51) is a consequence of plasma quasineutrality. 
It is useful to denote 


Lp = 


Zecj) 


( 2 . 6 ) 


where Z is the ratio of the ion to electron charge. 

Under these assumptions, we may integrate out the dependence of the ion distribution 


f The intricacies of the A:|| = 0 electron response (see section 12.1.11) do not affect the iner¬ 
tial-range theory to be presented here. 
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function on perpendicular velocities, so we introduce 

= — / d^v^Sfi, 

Til J 

and write the drift-kinetic equation for g in a 4D phase space: 
dq 

■^-hU||V||(g-H(/5FM) -Hu_l • V±g = C[g\ + x, 

[a 

ip = a 1 dv\\g, a = —, 

where -Fm is a ID Maxwellian with thermal speed Uth, 


FM = ^e ’'IIwth = 


0? 



uj_ is the E X B drift velocity, 


(5E X B - _ 

U_L = c —g 2 — = —^ ^ ^ 


(2.7) 

( 2 . 8 ) 
(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


C[g] is the collision operator and x ^ source term—both of which need a little further 
discussion, which we will provide in section r2.1.2l 

Note tha t while we will be referring to “slab” ion-temperature-g radient (ITG) turbu¬ 
lence (e.g., Cowley et al. 1991 : Ottaviani et al. 1997t Horton 19991) as the main physical 
instantiation that we have in mind of the kinetic problem described above, there will be 
nothing in our theory that would make it inapplicable to the ( i nertial range of) e lectron- 
temperature-gradient (ETC) turbulence ( Dorland et al. 20001 : Jenko et al. 2000ll . or in¬ 
deed to a generic case of electrostatic drift-kinetic turbulence with energy injection at 
long wavelengths. 


2.1.1. A nuance: Boltzmann closure and zonal flows 

In this context, we must come clean on an important detail. The Boltzmann closure 
(1^ for the electron density is, in fact, only valid for perturbations with /cy 0 be¬ 
cause it relies on electrons streaming quickly along the magnetic field lines to short out 



perturbations that have ky = fey = 0 (the “zonal flows”), namely, 


Sue 

Ue 


ejf - 4>) 

Te 


( 2 . 12 ) 


where f is the flux-surface average, which in our context is an average over y and z 
(formally, one gets this by first deriving the density response i n a toroidal, magnetically 
sheared system, as is done, e.g., in §J.2 of lAbel fc Cowlevll201,‘lL then taking the magnetic 
shear and curvature to be small and passing to the slab limit). 

This implies that, in order to find the (y, 2 :)-averaged (zonal) part of ip from the ion 
distribution function, at least the lowest-order finite-L armor-radius correctio n has to be 
kept (physically representing the polarisation drift; see iKromme 99.1 l 2 m nil . leading to 
a in equation dSlI) for the zonal part of ip having to be replaced with a = 2/kf.pf or 



(2.13) 
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In ITG turbulence far from marginal stability, these changes affect important quanti- 
tative details of the interaction between zonal flows and drift waves at the outer scale 
( Rogers et al. 2000t see also discussion around equation (12.431) 1. but do not matter for 
the inertial-range physics that we will focus on in this paper. 

A reader who is unconvinced may observe that equation (j2.9ll can be used without these 
modifications if, ins t ead of considerin g ITG turbulence, we consider ETG turbulence 
( Dorland et a/.|[2000l : ljenko et aLll2000ll . In this case, it is the ions that have a Boltzmann 
response (due to their large Larmor orbits, over which the density response from electron- 
scale fluctuations averages out), 

^ (2.14) 

Ti 

(which is = 5ne/ne by quasineutrality). The required modifications in equations 
and (1^ are 


ip ■ 


e(j) 

Ye’ 


a —>■ — 


ZTe’ 


Pi t Pe; ^th t Uthe — 


me 


(2.15) 


and ipFu —>■ —pFu in equation (12.81) . None of this affects anything essential in the 
upcoming theoretical developments. 


2.1.2. Injection, phase mixing, advection, dissipation 

The precise nature of the source term y in equation (12.8p will not matter in our 
theory, as long as it does not contain any sharp dependence on riy (i.e., is confined 
to low veloc ity moments). A random forcing is often a convenient choice for analytical 
theory ('e.g.. lPlnnkll2013 : Plunk fc Parkeiir2ni4 Kanekar et al. 2015 ). but a more physica l 
form in the context of electrostatic drift-kinetic turbulence in plasmas (e.g., HortonlllOOOl) 
arises from accounting for the presence of equilibrium density and temperature gradients, 
taken, conventionally, to be in the negative x direction^ 


X = - 


uj_ -ViniFu) 

n. 


PiVth dp 

\ 1 + 

f-l 

--V1 

2 dy 

Ln 


2 J Lt 


1 1 drii 1 1 dTi 

Ln Hi dx ’ Lt Ti dx 


(2.16) 


We shall see in section [27^ that these terms render the system linearly unstable and thus 
extract energy from the equilibrium gradients and inject it into the perturbed distribu¬ 
tion. 

The resulting perturbations are subject to two influences, linear and nonlinear, encoded 
by the second (uy Vyp) and fourth (uj_-Vj_p) terms on the left-hand side of equation (12.81) . 
respectively. The nonlinear term represents advection of the distribution function by 
the mean perpendicular flow, itself determined by the former. This involves coupling 
between different wave numbers and thus usually leads to spatial mixing (generation 
of small spatial scales) of the perturbed distribution. The linear term represents phase 
mixing—generation of small velocity-space scales in the perturbed distribution function. 
The simplest way to understand this is to notice that the homogeneous solution to the 
linear kinetic equation in Fourier space, dtg + iv^^k^^g = ..., is g ^ and the 

velocity gradient of that grows secularly with time, dy^^g = —ik\\tg. 

As fine structure in phase space is generated, there must be a means for removing it. 


t The erudite reader given pause by 1/2 rather than 3/2 in the prefactor of 1/Lt in equa¬ 
tion (|2.16ll will recall that we have integrated out the v± dependence. 
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This is why, even for a “collisionless” (meaning in fact weakly collisional) plasma, the 
collision operator C[g\ must be included in equation (I2.8|l . We hasten to acknowledge 
that, in pretending that the collision operator operates purely on we are ignoring 
that the vj^ dependence cannot in fact be integrated out of it; collisions will strive to 
isotropise the distribution and so the collision operator must necessarily couple vj^ and 
r;||. However, non-rigorously, when the collision frequency is small, 

ki\Vth, k±u±, (2-17) 


the collision operator’s essential contribution will be simply to iron out fine structure in 
velocity space and, given an initial distribution and a source that are smooth in v, only 
fine structure in uy ca n arise. Thus, it should suf hce to assume a simple model form for 
C[g]: for example, the Lenard fc Bernsteh] (1958) operator, 


C[9] = 


dv\\ 


IJL 

2 dv\\ 


+ vi\ g, U|| = 


n 

?^th' 


(2.18) 


The fact that this operator does not conserve mo mentum or energ y, while easily repaired 
if one strives for quantitatively precise energetics ( Kirkwo^ 194hl l. will not cause embar¬ 
rassment as collisions will only matter for high velocity moments (because large gradients 
with respect to ?;|| are necessary to offset the smallness of v). It is not hard to estimate 
the velocity-space scales at which collisions can become important: balancing C[g] ^ ujg, 
where w ^ k\\Vth and/or k±u± is the typical frequency scale of the collisionless dynamics, 
we find that the requisite velocity scale is 


5U|| 


(2.19) 


so the structure gets ever finer as —)► -fO. 

Finally, we are going to assume implicitly that equation (12.81) contains some regularising 
term to ensure a cutoff in k± —a necessity because of the spatial mixing associated with 
the nonlinear advection. Physically, the advection term will drive the system out of 
the domain of validity of the drift-kinetic approximation, to k±pi ^ 1 and larger. The 
precise way in which the energy is thermalised at these Larmor and sub-Larmor scales 
is a rich and interesting topic in its own right, involving a kinetic cascade in a 5D phase 
space (with nonlinear phase mixing in v \ now also occurring)—but these matters are 
outside the scope of this treatment fsee ISchekochihin et al. 2008 . 20091: Tatsuno et al 
2009t IPlunk et al. 2010l : Bahon Navarro et al. 20116). 


It is the competition between the two ways—linear phase mixing vs. nonlinear advection— 
of generating small-scale structure in phase space and thus enabling the energy injected 
by the source to be thermalised that will be the subject of this paper. 


2.2. Free energy 

We have referred to injection and thermalisation of energy many times now, and so 
defining precisely what we mean by “energy” has become overdue. 

Energy in Sf kinetics (i.e., in near-equilibrium kinetics) is the free energy associated 






























7 


Phase mixing vs. nonlinear advection in plasma turbulence 
with the perturbed distribution^ 


:F = -Y.^sSSs = Y.Ts6J d3v(/.ln/s)=E / 


( 2 . 20 ) 


where angle brackets denote spatial averaging and SSs is the mean additional (negative!) 
entropy associated with the perturbed distribution of species s. The last expression in 
equation (12.2011 was ob tained by letting fs = Tm .s +<5 f.s and expanding {fs In fs) to second 


order in Sfs (see, e.g.. ISchekochihin et aLll20(M note that {Sfs) = 0 because (fs) = Fms 


by definition). It is now not hard to establish that P = nipW, where 


W = 


dull 


2Fm 


2a 


( 2 . 21 ) 


is the quadratic quantity conserved by equation (12.81) . This can be shown either by using 
the Boltzmann-electron closure in equation (12.201) (viz., Sfe = {e(j)/Te)FMe, so the s = e 
term in gives rise to the ((p^) term in W) or directly starting from equation (12.81) . 
which gives us the following law of evolution of the free energy: 


dW 

dt 




(px) 



{gCjg]) 

Fm 


( 2 . 22 ) 


The first term on the right-hand side is the energy-injection term, which turns into the 
usual flux term for ITG (or ETG) turbulence if we substitute x from equation (I2.16P 
(see equation (12.371) below), and the second, negative dehnite, term is the collisional 
thermalisation of this energy flux. 


The Landau damping of the electrostatic perturbations is simply the transfer of free 
energy, via phase mixing, from the {(p'^) part of W to (5^)0: since ip = a J dv\\g, small- 
scale veloci ty-space structure in q is washed out in ip but of c ourse remains as free energy 
in (g^) (cl. iHammett fc Perkins 19901 : iHammett et Qi.lll99^ . 


t The understanding that this is the case can be traced back through a sequence 


of papers, from early, somewhat forgotten, insights to a more recent surge in ap 

precia- 

tion (iKruskal & ObermarJ 119581: iBernsteirJ 119581: iFowleil 119631. 119681: iKrommes & Hu 

ll994| 

Krommesll 19991: ISugama et a/.lll996l:lHallatschekll2004l:lHowes et a7ll2006l: ICandv & Waltzll2006l 

Schekochihin et a7l 120081. 20091: Scott! 120101: iBanon Navarro et a7ll20116llal: [Abel et al. 

boia 

Kunz et a7ll2015l: [Parker fc Dellai! 120151'). Note that we have not included in eouation (12.2011 


the energy of the electric and magnetic field, -|- (5B^)) /8tv (which is part of the general 

expression for the free energy; see, e.g., ISchekochihin et fflZjl200?j i because we are considering 
electrostatic perturbations (SB = 0) at scales much longer than the Debye length ((i?^) is 
negligible). 

j To be precise, from (p'^) and low-order (“fluid”) velocity moments of g to higher-order 
(“kinetic”) moments (see section [2.3.211 . 
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2.3. Hermite decomposition 

A natural way to separate the “fluid” part of the problem from the “kinetic” one and 
represent phase mixing is to expand the perturbed distribution in Hermite polynomials0 


5(^11) = - 9m, 


m—0 


V2^ml 


where Du = uy/uth and the “physicist’s” Hermite polynomials are 


iJ™(D||) = (-ire« 


dDn 


- Fuivu) = 


(2.23) 

(2.24) 

(2.25) 


The first three Hermite moments are the (ion) density (Sn), mean-parallel-velocity (uy) 
and parallel-temperature (5Ty) perturbations: 


ido(D||) = l 

Hl(Dy) = 2Dy 

^^2(Dy) = 4 fujj - 


5n if 

9o = — = —, 
n a 

= V2^, 

Vth 

1 ^Ty 

^ V2 T ■ 


9i 


Noting further that the source term, equation (I2.16|) . is 


X = 


PiVth dp 


Ho{v\\) H 2 {v\\) 


AiLj' 


Fu = 


Xo 


H2{v\\) 

2^2 


X2 


Fm 


(2.26) 

(2.27) 

(2.28) 

(2.29) 


2 dy Y Ln 

and that the streaming term in equation (12.8L uyVyg, couples Hermite moments of 
adjacent orders via the formula 

^ll^m(Dy) = ^ Hm+i(v\\) + mi7m_i(Dy), (2.30) 

we arrive at the following Hermite representation of equation (12.81) : 

—--futhVy—=yo = -.—p, —, 2.31) 

dt a uth 2L„ dy 

I + UX ■ Vx) a + „.,v„ (vi!s= + 2 ^) = (2.33) 

and, for m ^ 3, a universal equation retaining no traces of the temperature-gradient 
drive or Boltzmann-electron physics: 


d 


+ U_L ■ Vj_ 1 ffm + V y 


m -|- 1 


■ 9m+l 


gm-i = -vmgr, 


(2.34) 


at ' ■ "J 'MV 2 ' V 2 

^ Th is ha s attracted recur r iuR bursts of attention over many ye a rs, especial l y re¬ 
cently (iG radI [l9 4^: lArm stronel Il967l: iGra nt fc F e^ Il967l: lEltgrot hl Il974 ICrownfiddI Il977|: 


Harnmett et oZ.lll993l: iParker fc Caratill995l: lSmithlll997l: iN g et a7ll l999l: iWat anabe fc S^ugamal 


2004; Zocco_&_Schek^hlhln 20j4j;_ BlaHyejnj^l2oij;jLqu^^ et_alj^ 20b^ Hutch_et all l2013l . 


2014 : lKunk fc Parkeill2014l : Kanekar et aZ.ll201^ : lKanekarll2015l : iParker fc Dellarll2015fr 
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Note that we have taken advantage of the fact that Hermite polynomials are eigenfunc¬ 
tions of the Lenard-Bernstein operator ()2.18|1 . but ignored collisions in the m = 1 and 
m = 2 equations (this is allowed because we are assuming v —>■ -|-0 and so collisions will 
only be important at to ^ 1). 


2.3.1. Energy injection: slab ITG instability 

The first three equations are the standard three-field fluid system that describ es an 
ITG-unstable plasma at long wavelengths in an u nsheared slab (IC owlev et al\ 1991 ). The 
quickest way to obtain th e slab ITG instability ( Rudakov fc Sasrdeev 1961 : iGopoi et o,l. 


1967 : Gowlev et a,l\ 199lh is to balance the two terms on the left-hand side of equa¬ 


tion (I2.31|l . the first with the third term in equation (I2.32L and the first term on the left- 
hand side with the temperature-gradient term on the right-hand side of equation (12.331) . 
The resulting dispersion relation has three roots, of which one is unstable: 


a 





(fciiuth) 


2 / 3 , 1/3 




(2.35) 


where a;*T = kyPiVth/’^LT- This approximation is valid provided Ln/LT 1 and oj^t 
fc||Uth, although, as the growth rate grows with fey, the fastest growth is in fact achieved 
for fc||z;th ^ when the dispersion relation is a more complicated and somewhat 

unedifying equation. At fc||Uth the ITG mode is replaced by a sound wave, which, 

in a kinetic system, is heavily Landau dampedU 


2.3.2. Free-energy flows 

The temperature-gradient instability injects energy into the (/?, u\\/vth and per¬ 

turbations, all of which are comparable to each other in magnitude when fe||Uth ~ 
Because the three-field system is not closed^ there is a transfer of energy from SflflT to 
higher Hermite moments: the g^ term in equation (I2.33P provides the energy sink from the 
unstable (“forced”) moments and the g 2 term in equation (j2.34l) at to = 3 is the source 
for the higher moments; the energy thus received by them is eventually thermalised via 
collisions. 

To be more precise about these statements, let us rewrite the free energy (12.211) in 
terms of Hermite moments: 


Tj. 1 “k n / 2\ 

w = 


("ii) 


^th 


1 (ST^) 

4 7^2 


2 2^ 
m—3 


ial)- 


(2.36) 


t An elementary analy s is of t he slab ITG dispersion relation can be found, e.g., in Appendix 
B.2 of lSchekochihin et al\ ||2012| ~). Note that, in section [2.4.3l we will argue that the inertial-range 
fluctuations in fact have A:||Uth :^*t, and in sectional we will show that their Landau damping 
is suppressed in the nonlinear regime. 

J The only rigorous way to turn it into a closed system is to assume n ^ w, A:||Uth, k±u± in 
equation (12.3411 . whence gm g-m+i and so the heat flux is expressible in terms of the temper¬ 
ature gradient: VSga ~ (uth/v^z^)V||(;2 = (uth/2z^)V||(5T||/r. Putting this into equation (12.3311 
gives rise to a parallel heat conduction term. We are not, however, interested in this collisional 
limit. 
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Its “fluid” and “kinetic” parts satisfy® 


1 + a 




iV’^) 




‘"th 


4 T2 


2TLt 2T ^ II 11^^^’ 


= ^^|^('^r||V||g3) ->yf2 ^(Sm)- 


m—3 


(2.37) 

(2.38) 


m—3 


The first term on the right-hand side of equation (I2.37P is the injected energy flux and 
the second term on the right-hand side of equation (I2.38|) is the dissipation of that flux 
by collisions. In steady state, d(.. .)/dt = 0, we must have 

(^7||V||53) > 0 (2.39) 

because the collision term is negative-definite in equation (I2.38L and, therefore. 


{5T\\Ux) > 0 


(2.40) 


to achieve balance in equation (12.371) . The i nequality (12.391) implies a non-negative me an 
energy flux to higher Hermite moments fcf. iKrommes fc Hull 1994 : Nakata et al. 2012 ). 

How that flux is processed from being injected at m = 3 to being dissipated at m ^ 1 
(assuming i/ —> -|-0) is handled by equation (12.341) . This equation contains in a beautifully 
explicit form the two effects to which this paper is devoted: the phase mixing is manifest in 
that Qm is coupled to gm+i and Qm-i, providing a mechanism for pushing energy to higher 
m’s; simultaneously, all Hermite moments gm are advected (spatially mixed towards 
smaller scales) by the same fluctuating velocity uj_, determined, via equation (12.111) . by 
the zeroth Hermite moment, (p = ag^. 

2.4. “Fluid” turbulence theory 

Barnes et al\ ( 201ll) proposed a Kolmogorov-style theory of ITG turbulence, essentially 
ignoring the possibility of a leakage of free energy from the low Hermite moments to 
the high. W hile their theory is by no means unco ntroversial or the only offering on the 
market fe.g.. lGurcan et ol ]l20n9l:IPhmk' et al Jlmi, it does appear to match the results of 
numerical experiments (in the strongly unstable regime) and so it is worth both reviewing 
how it is constructed and examining to what extent it contradicts the statement made 
in the previous subsection that free energy mu st leak to higher Hermite moments. 

The scaling argument of Barnes et al. ( 20 111) addresses two main questions (as would 
any such argument aspiring to be a complete theory): 

(i) what is the effective outer (energy-containing) scale of the turbulence and the 
fluctuation level at that scale; 

(ii) what is the spatial structure of the turbulence in the “inertial range” between that 
outer scale and the small-scale cutoff? 


2.4.1. Outer scale 

The first question would be trivial for turbulence forced externally at some fixed scale, 
but for temperature-gradient-driven turbulence, relying on a linear instability, the system 
must decide where to have its energy-containing scale. This is not simply the peak scale 
of the growth rate, because at long wave lengths the growth rate will generally grow 
with wave number less quickly than will the nonlinear cascade rate (as we shall see in 

^ Note that both equations (12.3711 and (12.3811 will also contain sinks accounting for energy 
losses at small spatial scales. 
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section [2.4.211 and so in fact it is the largest scale at which there is 
w ill end up being the energy-containing scale. 

Barnes et al. ( 201ll) conjectured that this infrared cutoff will be 


parallel scale available to fluctuations: 


an instability that 
set by the largest 


fc||o 



(2.41) 


where Ly in our idealised homogeneous system is simply the parallel extent of the “box” — 
in a tokamak, it would be the magnetic connection length between unstable and stable 
parts of the plasma (Ly ^ qR, where q is safety factor and R the major radius). The 
perpendicular energy-containing scale is then given by 

— kyQPi— ~ /ciiQr^th ^ kyQPi ^ — (2.42) 

Lt L\\ 

because the instability would be supplanted by stable (in fact, Landau-damped) sound 
waves at smaller ky. Note that we require LT’/Ly <C 1 for the turbulence to occur in 
a scale range consistent with the drift-kinetic approximation k±pi <C 1. Finally, it is 
further conjectured that the zonal flows generated by the turbulence will have a typical 
she aring rate S'zf corn parable to the nonlinear decorrelation rate TjTi at the outer scale 
(cf. iRoeers et al.ll2000ll and, therefore, will isotropise the turbulence]^ 


kxo ~ SzpkyoTnio ^ kyo ^ k±o. 


(2.43) 


This is the only place in the theory where the zonal flows make an appearance as it is as¬ 
sumed that they do not completely dominate the nonl inear dynamics, in contrast to their 


alleged behaviour in the near-thre.shold regime ('e.g.. Dimits et al\ 20001: Diamond et al. 

2005 

2011: Giircan et a/.l 20091:iNakata et al. 2012: 

Ghim et aLll2013l: Gonnauehton et al. 

2014 

Makwana et al\ 2014). 


The energy-containing scale given by equation (j2.42|) , the amount of energy it contains 
is estimated by balancing the rate of injection by instability, a;*T, against the rate T~jg 
of nonlinear removal of this energy to smaller scales via advection by the turbulent flow: 


UJtT ~ k±opi 


Lx 


RlO 


k±ou±o ~ PiVthkxQtpo 


1 


Po 


kxoLr 


PtL\\ 

Ll 


(2.44) 


Finally, if one’s overriding practical concern is the calculation of the effective heat 
transport caused by the turbulence, one concludes from the above that the turbulent 
thermal diffusivity and the heat flux are 


-^turb 7 ^ Q ^ j-o • (2.45) 

/C_Lo ■L'T 

All of this is not particularly sensitive to the fact that, in making the argument that led 
to equation (12.441) . we completely ignored the possibility (in fact, the inevitability) that 


t It is possible to imagine (or conjecture) variants of drift-wave turbulence in which zonal flows 
are not strong enough to do this. In such systems, the saturated state at the outer scale is dom¬ 
inated by “streamers,” anisotropic st ructures with k^o - C kyp, whose radial extent is probably 
determined by the size of the system (Drake et aZ.lll98^ . Il99ll : ICowlev et oZ.lllfl^ : iRogers et al\ 
1 19981 : iDorland et oZ.|[2000l : iJenko et alll2000l l. In order for these structures to survive, they must 
be immune to the secondary ins tability that would otherwise give rise to zonal flows, which woul d 
in turn break up the streamers llRogers et oZ.l[2000 : lQuinn et aZ.l[2013l : [^nnaughton et a)]|2014fl . 
How a streamer-dominated outer-scale state channels its energy into an inertial-range cascade 
is not entirely well understood. However, we do not expect that the physics of this inertial range 
to be much different from that described below. 





























































12 


A. A. Schekochihin et al. 


some of the energy injected at the outer scale might be removed not by the nonlinear ad- 
vection, as if the system were purely fluid, but also by the phase mixing towards high m’s. 
The presence of such a transfer (of which there is , in fact, numerical evidence; see, e.g., 
Watanabe fc Sueama 20061 : iHatch et~d. 2OlT30; iNakata et al\ I 2 OI 2 II would only break 
our argument if the rate ^ fciio'i^th of this transfer were substantially larger than the non¬ 
linear advection rate and so if the dominant balance were ~ feyo'^^th ^ k±QU±o. But 
this is obviously impossible as one cannot saturate a linear instability by a linear mecha¬ 
nism: there would not be anything in the theory to determine the saturated amplitude^ 
In view of equation (I2.42L the phase mixing rate is, in fact, of the same order as both 
and k±ou±o. Therefore, it cannot affect the basic scalings—although for the purposes of 
quantitative transport modelling, it is quite crucial to know by what fraction of order 
unity it might cut the nonlinear mixing rate, a key preoccupation in the development of 
“Landau-fluid” closures for plasma turbulence in fusion contexts faammett et al. B)92 , 
I 993 I : Borland &: Hammett 1993 : Beer fc Hammett 1996t Snvder &: HammettI 2001^ 


A question that is much more sensitive to whether phase mixing is nonnegligible is the 
structure of the inertial range. 


2.4.2. Inertial range: perpendicular spectrum 

How is the energy inject ed at (fcj_o,^||o) cascaded to smaller scales? Ignoring phase 
mixing, Barnes et al. ( 201 il l proposed to calculate the dependence of the turbulent am¬ 
plitudes on scale via the Kolmogorov assumption of constant energy flux: at each scale 
energy ipA is transferred (locally) to the next smaller scale over the cascade time 


fl 

T-nl 


k±u±<p^ oc k\ip^ = const 


7 -2/3 
ipock^ , 


(2.46) 


where we used u±_ oc k_\_I> (see equation ()2.11|l b Note that, both here and in similar 
arguments that will follow, we do not make a distinction between the energy content of 
low-TO moments, assuming 

u\\ STu 

— (2.47) 

and possibly also ^ a few more low-m moments of g, although we do assume that there 
is not a substantial energy leakage to asymptotically large m’s. The ID (perpendicular) 
spectrum is then 

(2.48) 


/ 2 

dfc||(|v?kP) ^ 


-7/3 

_L ’ 


where (• • •) now denotes a time or ens e mble average. This scaling is supported both by 
numeri cal simulations of iBarnes et al. d^llll and , apparently, by those done by other 
groups ( Hatch et a7ll2013Ll20l4 Plunk et aL[2015L who confirm finding the same scaling, 
without, however, providing plots). 


2.4.3. Critical balance 

The structure of the turbulence in the parallel direction can now be inferred via a 
causality argument known in th e astrophysical MHD literature as “critical balance” 

( Goldreich &: Sridhar 1995 , 1997 ; BoldvreT^ 2005lj and e merging as a universal scaling _ 

principle for strong turbulence in wave-supporting systems ( Cho fc Lazarianl2004 : Schekochihin et al 


t Again, focusing on turbulence far above the threshold, we are going to ignore the possibility 
of a more sophisticated scheme involving zonal flows. 

t This is because the typical rate for coupling these moments is fc||Uth, which will shortly be 
argued to be comparable to the nonlinear rate at which these moments change, equation (12.501) . 
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2 OO 9 I: Nazarenko fc Schekochihin 2011 k fluctuations cannot stay correlated at parallel 


distances longer than those over which linear communication happens at the same rate 
as the nonlinear decorrelation: thus, fluctuations are uncorrelated for 

4/3 


^11 ^th "S k±^U±^ oc 


k\\L\\ 5 , 


f—V 


(2.49) 


Here and in what follows, we shall adopt a nondimensionalisation Ly = 1 and fcj_o = 1, 
so the above condition will henceforth be written fey < 

This argument implies that, at any given fc_L, the “energy-containing” parallel scale 
will be given by the “critical-balance” wave number: 




fell 


l4/3 


(2.50) 


another scaling that was confirmed numerically by iBarnes et al\ (|2011r) . The consequent 
scaling of the ID parallel spectrum is (using equation (j2.50l) in equation (j2.46l) fi1 


oc fcy ^ =2tt j Ak±_kx{W-k\^) 


k\\ 


oc fell 


(2.51) 


Note that, under this scheme, the drift waves are slow in the inertial range because 
oc ky whereas k±^ux oc so the relevant frequency in equation (12.491) is indeed 
~ fcyuth, not uJtfT- By the same token, energy injection by the temperature-gradient 
instability is slow compared to the nonlinear cascade rate, so, effectively, the instability 
only operates at the outer scale, while the fluctuations that carry the injected energy 
through the inertial range are more akin to ion sound waves than to drift wavesQ 


2.4.4. Constant flux is inconsistent with robust phase mixing 

In section 021 we will explain how to derive from these arguments the scaling of the 
2D spectra for any fc_L and fcy. However, we must first discuss the key point that the 
constant-flux assumption (12.461) cannot be consistent with both the idea that the energy 
resides along the “critical-balance curve” (12.501) and with phase mixing taking energy out 
to large m’s at the rate ^ fcyUth—simply because the latter would mean that the energy 
in the low-m moments is not conserved and so need not be fully transferred nonlinearly 
to smaller scales. 

The simplest way to explain the implications of this for the spectra is to replace 
equation (12.461) by a simple mock-up of an evolution equation for E^(k±) (cf. iBatchelor 
1953HHowes et aLll20n8^ : 


dE^ 

dt 


de 

dk± 


- lE^ 


kxE^ 

T„1 


7 ~ ^ f 


-1 

nl 


kl 


\J kxE^, 


(2.52) 


f Another way of arriving at this spectrum and at the critical balance (iBeresnvald 120151 ') is 
to start with the constant-flux conjecture applied to the scaling of amplitudes with frequencies, 
rather than wavenumbers: gfoj ^ const ip oc (ICorrsirJll963f) . The frequencies of the 

perturbations will be ta ~ fciiri^h . hence the p a rallel scaling (12.5111 . 

j This also explains why the Bames_et_ad (DOllh cascade should asymptotically override the 
nonlinear transfer proposed by Giircan et aZ.I ([20091 ): the latter authors argue, effectively, that 
the cascading of the energy to small scales is done by the nonlocal shearing of the drift waves by 
zonal flows, which they assume to occur at the rate ~ fcxMzF, where uzf is a scale-independent 
zonal velocity; this, via a constant-flux argument analogous to (12.461) . gives oc kj^. However, 
if the zonal shearing rate is comparable to the energy-injection rate at the outer scale (which 
we also assume; see equations (|2.43ll and (I2.44|l : k±oUzF ~ Szf ~ ~ w*t), then it will 

be smaller than k±u± for fcx > fcxo. Note also that a nonlinear transfer rate oc k± could not 
effectively dominate the injection rate, lo^t, which is also oc k±. 
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where e is the energy flux and 7 is the effective rate of phase mixing (Landau damping), 
which, by the critical-balance conjecture ()2.50|1 . is of the same order as the cascade 
rate . Assuming steady state in equation (12.521) and letting yrni = C = const ^ 1 
(independent of fcj_, as per critical balance), we get 

de ^ 




£ OC 


E^{ks_) OC k 


-(7+2C)/3 


(2.53) 


dk± k± 

Thus, the flux decreases with increasing wave number and so the spectrum is steeper 
than the constant-flux solution (j2.48|) . The power laws that emer ge in such dissipative 
systems are generally hard to predict and probably nonuniyersal ('cf iBratanoy eA, Qi.ll2013l : 
Passot fc Sulemll2015ri —in our case, because they depend on an order-unity prefactor (^) 
in the critical-balance relation (12.501) . rather than on some dimensionally and physically 
inevitable scalinglj However, numerical—or, indeed, experimental—evidence does not 
appear to support spectra th at are significantly steeper than k~^ at long (above the Lar- 


mor s cale) wavelengths (e . g.. iHenn eanin el oi.ll2004tKtorler <V: ,lenlcollZl)U»t itJasati et al. 
2OO9I : Vermare et al. 20 lit iBarnes et Qi.ll201ir Kobavashi fc Gii rcanll201 5ll. Furthermore . 


y steeper tnan fc , at long (above tne l^ar- 
j] 2 OO 4 I: Gorier fc Jenkcj 20(18 : Casati et al 


numerical investigations by Teaca et al\ ~i 201^ 2014[l and Bahon Navarro et al. ( 2014 ) 

confirm local nonlinear energy transfer and possibly even constant fluxes, albeit with 
a number of caveats regarding non-asymptoticity of the simulations, consequent possi¬ 
ble non-universality of their results, as well as distinctly measurable, if not dominant, 
amounts of dissipation (meaning, in their context, phase mixing) everywhere. 

In what follows, we shall see that, in a sufficiently collisionless plasma, the constant-flux 
assumption is safer than it might appear. 

2.5. Hermite “cascade” 

As the last bit of essential background, let us consider what happens with free energy in 
phase space if we treat phase mixing as the dominant process and ignore nonlinearity— 
the opposite extreme to that pursued in section [2j4l 

Returning to equation (j2.34j) and dropping the advection term uj_ • Vj_ for the time 
being, we perform a F ourier transform in the paralle l direction and introduce the following 
very useful functions ( Zocco fc Schekochihinll201ll) : 


ffm(fc||) = (isgnA;||)™ 5 „^(A;||), (2.54) 

where 5 m(fc||) are the Fourier-Hermite harmonics. The (linearised) equation (12.341) then 
becomes 

dgm , 1^11 kth 


dt 


V2 


'-{Vm+ Igm+i - y/rngm-i) = -vmgr, 


(2.55) 


f It is easy to see that ^ < 1. Indeed, the nonlinear cascade rate that follows from equa¬ 
tion (12.531) is k±u± OC , which can only overcome the injection rate associated with the 

temperature gradient if ^ < 1 (see discussion at the end of section 12.4.311 . The extreme case 
^ = 1 gives oc kj^. One can obtain such a spectrum if one assumes that the fluctuation 
energy present at each scale, not just at the outer scale, is determined by the balance between 
the instability growth rate, the nonlinear decorrelation rate—and also the phase mixing, which 
removes the energy to high m’s, so there is no need for a constant flux. Then each scale be¬ 
haves as the outer scale described in section [2.4.1l ('co oc as in equation (12.4411 ). We consider 
this scenario much too fanciful (it would require quite a complicated set of arrangements in 
the (fcj_,fc||) space) and rather unlikely for a system far from the threshold. Note also that the 
restriction ^ < 1 would not apply in a system where the energy injection rate is not proportional 
to k±, e.g., one where y in equation (12.81) is just a large-scale force and so the injection occurs 
only at the scale of the force. Then the non-universal spectrum (12.5311 can be steeper than kj^, 
although we must have ^ < 4 in order for k±u± to increase with fcx and so for the nonlinear 
transfer to stay local. The steepest possible spectrum in this case is, therefore, oc kj^. 
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The point of these manipulations is that they have made the phase-mixing term on 
the left-hand side of equation (j2.55|) look like a derivative with respect to m. Indeed, 
assuming, in the limit of to ^ 1 , that we can treat pm as though it were continuous and 
differentiable in to (an assumption that will come under close scrutiny in section IQ) , 
i.e., gm±i ~ 3m ± dmgm, we have 


\/m + Ig-m+i - Vmgm-i = 



1 

1 T 3m-t-l 9m—1 
TO 


j 9gm 
dm 


dm ^ 


Thus, equation (|2.55ll becomes 


dt 


V2 |fc||| 


vthm 


1/4. 


dm 


rn^^^gm = 


-vmgr, 


(2.56) 


(2.57) 


Introducing the Fourier-Hermitefree-energy spectrum Cm(fc||) = {\9m{k\\)\^) = (|3m(fc||)P), 

we find 

—^ +Cm =-2i^mCm- (2.58) 

In st eady state, the solution is ( Zocco fc Schekochihin 2011 : cf. Watanabe fc Sueamal 
2004 

VTO \ 2 v 2 i^ / 

where 4(fc||) is the constant of integration. Below the collisional cutoff, to -C TOc, the 
power-law scaling Cm oc m~^l'^ is the solution corresponding to constant free-energy flux 
in Hermite space (the Hermite flux is the expression under dm in equation (I2.58|) b It 
is possible to show quite rigorously (by direct Hermite transformation of the Landau 
response function) that this is indeed the Hermite-space solu tion that arises in a linear 
system with external forcing at low to and Landau damping ( Kanekar et a7^20I5^ . 

The solution (12.591) has two important properties. Firstly, the free-energy dissipation 
associated with it (the last term in equation (12.381) ') is dominated by Hermite moments 
with TO ^ TOc and does not explicitly depend on the collision frequency (assuming A(fc||) 
does not), 


D = v E W(3m) = V EE mCm{k\\) ~ u E/ <irnmCmikii) = Y^ 


\k\\\vthA{kii) 

V2 ■ 


(2.60) 

Secondly, the total amount of free energy stored in the phase space in order to achieve this 
finite dissipation (corresponding to finite amount of injected power) diverges as z/ —0: 


W 


E / d"lC'm(fc||) = ^ 


r(i/3) 4(fc||) 

3^A\/2 (|fc|||uth)^/^z^^'^^ 


(2.61) 


( Kanekar et al ][20il). 

Thus, if we thought that Landau damping in a turbulent system works in the same 
way as it does in a linear one, we might have to conclude that, rather than staying in low 
m’s and being nonlinearly cascaded to small spatial scales, as in a fluid problem, the free 
energy fills up phase space and diss ipates on coll i sions. A ded icated study of the Hermite 
spectra of slab ITG turbulence bv iHatch et a,l\ ( 20131 [20l4 showed that this does not 
happen, with Hermite spectrum of the free energy following a much steeper power law 
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than equation (I2.59|) and the wavenumber spectrum consistent with equation (I2.48[) . In 
what follows, we will show how such a solution can emerge Isection [4.4.21 has the answer 
and appendix E] the physical basis for it; see section 14.71 for the noirlinear versions of 
equations (12.611) and (12.601) 1. 


3. Formalism 


3.1. Phase mixing and anti-phase-mixing 


Our first order of business in constructing an appropriate irrathematical description for 
phase-space turbulence is to reexamine our rather blithe assumption in section 12.51 that 
gm{k\\), defiired by equation (12.541) and satisfyiirg equation (12.551) (to which the nonlin¬ 
earity will be restored in section IQ) . cair be treated as continuous in m. 

Consider 


1 < TO <C 



2 


(3.1) 


If we assume that the rate of change of gm is small compared to v^l^ykth, equa¬ 
tion (I2.55D tells us that, to lowest approximation, 


\/m + Igm+i - Vmgm-i = 0 


~ Qm—1’ 


(3.2) 


This has two solutions: 


so, in fact, either g^ or {—l)^g„ 


ss zLg 

mi (3.3) 

can be treated as continuous in to. We therefore 


2015 


intr oduce the following d ecomposition (which we already used in iKanekar et al. 
and Parker fc Dellar 2015 ) 

+ (3'4) 

where 

9m T 9m-i-l - / \rn9m 9m-i-l 


9m 


9m = (-1)^ 


(3.5) 


2 ’ V / 2 

can both be assumed continuous in to. Evolution equations for these two types of modes 
can be derived by adding or subtracting evolution equations (12.551) for g^ and gm-i-i and 
then expanding in large to in the same fashion as we did in section [2.51 The result is 


^ ± |fc|| |uth TOl/4^TOl/4g± = 


-vmg„ 


(3.6) 


Manifestly, the “-I-” modes are the phase-mixing modes, propagating from small to large 
TO, whereas the ” modes propagate from large to small to and thus represent “anti- 
phase-mixing” : free energy co ming back from phase s pace , a poss i bility earlier mooted, in 


pnase-mixmg : tree energy co ming bacK: trom pnase s pace , a poss i Diiity earner mootea, in 
somewhat different terms, bv iHammett et a)T ( 19931) and Smith ( 1997 ). We shall discuss 
the energetics of this process more quantitatively in section 151^ 

In a linear problem, in the absence of free-energy sources at high to, the only solution 
that satisfies the boundary condition gm^oo —t 0 is 5 “ = 0 , so there will be no anti¬ 
phase-mixing and the treatment in section 12.51 is correctly As we are about to see, the 
situation changes once nonlinearity is accounted for. 


f IKanekar et all (I2015D showed that in a (forced) linear problem, the spectrum of the ” 
modes is oc and so subdominant to the spectrum (12.591) of the “-I-” modes. This does not 

mean that there is some small subdominant amount of anti-phase-mixing in a linear system, but 
is rather due to the interpretation of and as being forward and backward propagating 
modes in m space being correct only to lowest order in 1/m. Note that this interpretation breaks 
down also at such large m that the inequality dni is no longer satisfied. When m (|fc||nth|/i^)^, 
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3.2. Nonlinear coupling and plasma echo 

Let us now restore the nonlinear advection (the second term on the left-hand side of 
equation (12.341) ') and Fourier transform it in the parallel direction: 

= - [uj-■(fc||) =-^ uj_(p||) • V_Lg™(g||). (3.7) 

^ P||+9||=fc|| 

Then the nonlinear term that must be added to the right-hand side of equation (|2.55p is 


dgm 


nl 


= -(isgnfc||)™ [u_L • V±gra] (fc||) = U_l(P||) • V_L 

P||+9n=fc|i 


/ sgnfcii Y 

Vsgn9|| J 


9m{q\\)- 


(3.8) 

Finally, adding or subtracting the above for the m-th and (m -I- l)-st Hermite moments, 
and using the decomposition (13.41) . we find the nonlinear term for equation (13.61) : 


dt 


± ^2 |fc|| \vth + umg^ = 


dm 

'^j-(7’ii) ■ 

P||-l-9||=fe|| 


St „Jtiiq\\) + Sk 


.9ir 


.911 


ffm(9||) 


(3.9) 


where = [l ± sgn(fc||g||)] /2, i.e., is non-zero (and equals unity) only if /cy and gy 

have the same sign, and S~ is non-zero (and equals unity) only if they have the opposite 
sign. 

The key development manifest in equation (13.91) is that the advecting velocity field can 
couple parallel wave numbers of opposite signs and thus produce anti-phase-mixing ” 
modes out of phase-mixing “-I-” ones and vice versa; = 0 is no longer a valid solution. 
This is a manifestation of the textbook plasm a-physics phenomenon known as plasma 
echo ( Gould et al. 1967 : Malmberg et a/.l[l968ll . The importance of it in our discussion is 
that once the free-energy flux through phase space is not compelled to be unidirectional 
towards high m’s (as it was in the naive treatment of section [23]) . all bets are off as 
to the effectiveness of Landau damping/phase mixing as a dissipation mechanism in a 
nonlinear system. 

3.3. Dual kinetic equation in phase space 

Equation (13.91) can be recast in a remarkably simple form if we introduce a change of 
variables and a rescaling of git'- 


/(s.^ll) = 


-^1/4, 


gtiih) 


if 

if 


fcy ^ 0, 
fcii < 0. 


(3.10) 


For any given fcy ^ 0, the original distribution function is reconstructed in the following 
way, via equations (12.541) and (13.41) : 


gm(fcy) = (-*)™ hVm, fcy) -f (-1 )”"/(Vto, -fc||) 
The new function / satisfies 


ffm(-fc||) =5m(fc||). (3.11) 


df hvti^df 
dt ds 


p\\ 


V_L/(fcy -p\\). 


(3.12) 


the collisional term in the right-hand side of equation (12.551) is dominant and the solntion 
is grn ~ {\k\\\vth/v\/2m)g.m-\ <C gm-i- Therefore, in this approximation, the two modes are 
gti ~ ffwi/2 « (—l)"‘gi7i and so they formally have the same energy. 
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The echo effect in this equation looks explicitly like mode coupling from positive to nega¬ 
tive parallel wave numbers, or vice versa, whereas the phase mixing and anti-phase-mixing 
are simply propagation in s with velocity k\\Vih_/\/2. We will make repeated references to 
this equation in the scaling arguments of section 01 

Equation (|3.12l) is a kinetic equation in phase space dual to the original kinetic equation 
(|2.8I1 . with the variable s (or \/2s/vth) effectively acting as a Fourier dual to riy—this is 
not a huge surprise because for m ^ 1, Hermite polynomials are well approximated by 
trigonometric functions in r;||, with “frequency” •\/2m/rith: 

Hm{vi\)e RiV2(—j cosi—— r;||- —j. (3.13) 

It is worth stressing that, while the functions &re subject to reality conditions, 

inherited from gm via gm (see definitions (I2.54|) and (I3.5|) l. 

9mi-ki\) = g^{ki\) ^ gm{-ki\) = g^{kii) g^{-k\\) = [gt,{k\\)]* , (3.14) 


the function /(fey) has no such property because it has been spliced together from the 
positive-fc|| values of g^ and the negative-fc|| values of 5“ and there is, a priori, no 
symmetry between the “-I-” and ” modes. 

Let us reinforce this point by showing that a solution of equation (13.1211 can only have 
the property 

/(-fell) = r (fell) ( 3 . 15 ) 

if the phase mixing is ignorable (this is worth noting because if equation (|3.15p does 
hold, then the free-energy flux in Hermite space vanishes, as per equation p.27|) : we 
will make good use of this argument in section Ol) . Taking the complex conjugate of 
equation (13.121) and subtracting from it the same equation written for /(— fey), we get 




r(fcy)-/(-fcy) 


k\\vth d 


r{k\\)+fi-k\\) 


pii 


ds . 

r{k\\ +P\\) - /(-fcy -py) 


(3.16) 


where have used u^(py) = uj_(— py) and then changed the summation variable py —>• —py 
in the sum involving f*. Equation (I3.16|) is compatible with the condition (13.1511 only if 
the phase mixing term can be ignored—which might happen because fcy is small and/or 
because / depends on s in such a way that the phase-mixing term is subdominant at, 
say, high s. 


3.4. Free-energy spectrum and free-energy flux 

Since we are going to discuss free-energy spectra and free-energy fluxes in phase space, 
let us provide the formal definitions and evolution equations for them. 

We define, in the same way as we did in section [231 


C™(fcy) = (|p™(fcy)P) = (|p™(fcy)P). (3.17) 

Then, using equation (I2.55|) with the nonlinear term given by equation (13.81) . we have 

^ + + = (3.18) 
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where the free-energy flux in Hermite space is (cf. IWatanabe fc Sugam ^ I2nn4 

rm(fc||) = v^2(to+ 1) |fc|| |'(;thRe(gm+i5^) = A/2(m"+T) ki\VthlTn{gm+igtn) ■ (3-19) 

These expressions are exact. In the limit of to ^ 1, both Cm and Fm are continuous in to 
( even if gm alternates sign, equation (13.31) 1. so we may rewrite equation (13.181) as follows 


dC„ 


dTr, 


+ 2vmCm = 


dC„ 


nl 


dt dm V dt 

Using the definition (13.51) of the “±” modes and defining their spectra 

we notice that, still exactly, for any to, 

Tm = \/2(to+ 1) Ifeyluth (C'+ - C^) . 


(3.20) 


(3.21) 


(3.22) 


Thus, the Hermite flux is exactly proportional to the difference between the spectra of the 
“+” and ” modes. The sum of these spectra is the free energy, but only approximately, 
for TO 1: 

Cm T Cm+1 


Ct + C- = 


c„ 


(3.23) 


The evolution equations for Ci^ at large to are best obtained from equation (13.91) : 
-^±1^11 kth + 2i'mC^ = 


-2ReJ2 ([5m(fc||)]*Uj.(p||) • V_L ['5^|.,„5m(9||) +'5fc||.9„5m(9||) 

P||+9||=fc|| 


(3.24) 


The sum of these two equations gives us back equation (13.201) with Tm given by equa¬ 
tion (13.221) (with TO >> 1). Another, more compact, way to write equation (13.241) is in 
terms of the spectrum of the function / introduced in section 13.31 Defining 

F(s,fc||) = (|/(s,fc||)|2), (3.25) 

we infer from equation (13.121) : 

% + (.r(fc||)u^(p||). V^/(fc|| -pii)) . (3.26) 

pii 

Note that, whereas C'm(fc||), C^{k\\) and rm(fc||) must all be even in k\\ because of the 
reality conditions (I3.14p . there is no such constraint on F{k^^) and, in fact, it is the odd 
part of A(fc||) that sets the Hermite flux: in view of equation (I3.22p . 


r^(/c||) Ri v^|fc|||uth [A(s,|fc|||) - T(s,-|fc|||)] . 


(3.27) 


The next step in the formal solution of the problem is to solve equation (13.261) for 
7^(5, fc||). However, even in principle, this is only possible if a suitable closure is found 
for the tri ple correlator in th e right- hand side. A particular solvable model will be dis- 


lor tne tri ple correlator m tn e rignt- nanQ siOe. A particular solvable model will be dis¬ 
cussed in ISchekochihin et ( 20161) . but it will come at the price of decoupling the 


advecting velocity from the advected distribution function (i.e., considering a “kinetic 
passive scalar”, rather than the fully self-consistent turbulence problem). In general, as 
always with turbulence problems, we are reduced to (or blessed with) having to resort 
to phenomenological scaling theories, which we will pursue in the next section. 
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4. Scaling theory 

In constructing the scaling theory for our turbulence in phase space, we shall continue 
to consider as sensible and valid the arguments in section [2.4.11 that led to estimates of 
the outer scale (equation (I2.42|) l and the amplitude of ip at that scale (equation (12.441) '). 
We will, therefore, focus on what happens in the inertial range. 

The argument that is presented in this section is quite long because, even within the 
inertial range, the phase space splits into several regions, where different physics are 
at work—and building the full picture involves investigating each of these regions and 
matching free-energy spectra at their boundaries. A road map to what is done where is 
provided by the subsection headings and by the overall summary in section 15.11 which 
an impatient reader might find it useful to read first. 

In what follows, wherever our expressions appear to be dimensionally incorrect, this is 
because the wave numbers are normalised to the outer scale: 

^ = fc||L|| -)>/cy, -t fc_L; (4.1) 

'C||0 ^-LO 

we will also omit, wherever this makes exposition more rather than less transparent, 
such dimensional factors as Uth, Pi, etc. We remind the reader that at the outer scale, the 
parallel-propagation/phase-mixing and the nonlinear-advection time scales are assumed 
comparable, feyoUth ^ 


4.1. Spectra in the phase-mixing-dominated region 

The first two terms on the left-hand side of equation (13.121) describe propagation of a 
perturbation in s with time, along the characteristic 


s = 


k\\Vth ^ 


(4.2) 


If we consider fcy > 0, perturbations will phase-mix in an unfettered way for at least a 
time comparable to the time it takes the nonlinearity to couple these perturbations to 
different wave numbers: 

ik±u±)~^ oc kj^, (4.3) 

where r is the scaling exponent of the nonlinear decorrelation rate. This means that 
whatever spectrum, denoted Ey,{k\\ ,k±), prevails at low s (and so low m), it will simply 
propagate to higher s as long as 


S < k^VthTnl 


k\\vth 

k±u± 



(4.4) 


We can rearrange this statement to mean that, for any given m = the part of the 
wave-number space satisfying 


k\\Vth ^ ^/mk±u± ^ fcy > ^/mk'^ 

will contain an exact replica of the low-m spectrum^ 


Ej{s, |fcy|,fc_L) = Vm£’+(fcy,fc_L) - E^{kj_,k\\) ~ k’]_k^^ “. 


(4.5) 

(4.6) 


t We assume that there is no discontinuity in the Hermite spectrum at low m, i.e., that the 
low-s limit of the solution to equation (13.1211 (which is technically only valid for s = ^/m ^ 1) 
will smoothly connect onto the spectra of low-m “fluid” moments ip, uii, ^Ty, etc. and also that 
the spectra of these quantities all have the same scaling with fcy and k±. 
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Here we have defined 2D spectra 

Ef{s,k\\,k±) = 27rfc_L(|/(s,/c||,k_L)P), 

S™(fc||,fc_L) = 27r/cj_(|g±(fc||,k_L)P), (4.7) 

E^{kii,k±) = 27rfc_L(|¥5(fc||,k_L)p), 

where (/7(fc|| ,k_L) etc. are Fourier transforms of the original fields in all three spatial 
directions. We shall refer to the lower bound on fey (or upper bound on k±) defined by 
the condition (14.5L fey ^ y/mk"^, as the phase-mixing threshold. 

The scaling exponents a and b in equation (14.61) are as yet unknown. One of them, b, 
can be determined in a purely “kinematic” way: since it describes the low-fcj_ (see equa¬ 
tion (14.51) 1 asymptotic behaviour of the spectrum, it must, in a homogeneous isotropic 
system, be 6 = 3 (the derivation of this result, which is quite standard, is given in ap¬ 
pendix |X]— it describes the spectrum at perpendicular wavelengths that are longer than 
the perpendicular correlation scale of perturbations with a given fcy). 

Thus, we have found a phase-mixing-dominated region (as we shall henceforth call it) 
of the phase space, with spectra 

^3 f^-a 

■£^m(^lh^-L)- E-{kll,k±) E+{kll,k±), kll>y/^kl. (4.8) 

v 


These and all subsequent spectra that will emerge are sketched in figure O which the 
reader is invited to consult for illustration (and preview) of the upcoming results, as they 
emerge. 

Unsurprisingly, in equation (14.81) we have a Ijy/m Hermite spectrum—the standard 
linear result already derived in section [5^ The anti-phase-mixing component of the free 
energy (U“) must be small compared to the phase-mixing one here because this is the 
part of phase space where the nonlinearity has no time to exert any influence and so 
there will not be any echo effect. 

While we do not yet know the exponent a (it will be deduced, in two different ways, 
in sections [4.21 and l4^ . it is clear that U+(fey,fc_L) must decay sufficiently fast with 
fcy in order for the total free energy not to diverge at short parallel wave lengths. 
This tendency for the free-energy spectrum to decay sharply at parallel wave numbers 
bounded from below (or, equivalently, at perpendicular wave numbers bounded from 
above) by the phase-mixing threshold fcyuth ^ k±u± (one might also call this thres hold 
the “phase-space critical ba lance”) was recently reported bv iHatch et all ( 201311201^ (cf. 
IWatanabe fc Sugamal 2004t) in their simulations of slab ITG turbulence (they, however, 
had a different explanation for it). 


4.2. Spectra of low moments 

As we explained in section 14.11 the spectrum (j4.8|l is inherited (propagated by phase 
mixing) from low m’s, so we must have 

E^{ki\,k^)^klkf, fc||>fcl. (4.9) 

Thus, this is the 2D spectrum of the electrostatic turbulence on the short-parallel- 
wavelength side of the critical-balance condition (j2.50ll . 

As we argued in section [2.4.31 the critical balance is essentially a causality condition and 
so the spectrum at the long-parallel-wavelength side of the critical balance, /c|| < k'^, must 
reflect the fact that the perturbations at these parallel scales are essentially uncorrelated. 
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The spectrum of such uncorrelated perturbations is the spectrum of white noise, so 

fc||<fcl. (4.10) 

Matching this with equation (14.91) along the curve fey ~ k^j_ gives 

a=^. (4.11) 

r 

If a > 1, then k\\ ~ k^^ is the energy-containing parallel scale for any given k±. The 
ID perpendicular spectrum is, therefore. 


E^{k^) = / dfc|| £:,^(fc||,fc_L) 


d/CM fcn/c 


0;„-c 




This immediately implies a consistency relation between c and rl3 

~ ~ k^ip ~ k']_{k±E^Y^'^ => r = 5 — c. 

Finally, the ID parallel spectrum for any given ku is dominated by k± ^ kJ : 


El{kY = J dfci £;^(fc||, fc±) - dk^ k^Yl" ^ k 


-(c-l)/r 


(4.12) 

(4.13) 

(4.14) 


4.2.1. Scaling exponents under constant-flux conjecture 

Note that so far, we have invoked no cascade physics, but, in order to determine the 
exponent c, we do now need to make an assumption as to how energy is passed from 
scale to scale. Energetically, only the wave-number region fey < k'^ matters because at 
larger fcy , the spectrum is assumed (and will be confirmed) to have a steep decay with fcy 
(equation (14.91) 1. We shall call it the advection-dominated region and anticipate that phase 
mixing there will not be a significant energy sink, i.e., the anti-phase-mixing energy flux 
due to the echo effect will on average cancel the phase-mixing flux, leading to effective 
conservation of (</5^). Then we can return to the constant-flux argument of section [2.4.21 


~ const k^EY{k_i) ^ ^ c = 1-I-2?’, (4-15) 


where equation (I4.12|) was used to obtain the last relation. Combining equations (14.151) . 
(14.131) and (14.111) . we find 


11 


3’ 


0 = 5. 


(4.16) 


This gives us back the lBarnes et al\ (j201l[ ) ID spectra: 


EY{k^)^kfY/\ ElikY^k-^ (4.17) 

via equations (14.121) and (I4.14p , respectively. We have now also learned what the full 2D 
spectrum behind these ID ones is: combining equations (14.91) and (14.101) with the scaling 


t We remind the reader that p> here and in all similar calculations in this paper is not the 
Fourier transform of the potential, but rather its amplitude corresponding to the scale k~Y (this 
can be thought of, for example, as the typical magnitude of the potential’s increment across a 
distance ^7^)- relationship to the Fourier transform ip-k. and to the ID spectrum Ef{k±) was 
given in equation (12.481) . This can be understood dimensionally or by noticing that the energy 
associated with a given scale is the integral over the energies contained in the wave number 
k±_ and larger, ~ dkfE;!;(kf) ~ k±E;^{k±) (the latter relation holds as long as the ID 
spectrum is steeper than A:7^). 
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(a) 


E4,ik^\,kPj 



E^{k\\,kj_) 



(c) 


Figure 1. Spectra of low Hermite moments, Etp{k\\,k±)-. (a) in the (fc^,fe||) plane, (c) vs. fc|| at 
constant k±, (b) vs. k±_ at constant fey. All plots are logarithmic. 


exponents (14.161) . we have 

Eip(k\\ , kx) 


if 

k\k-^ if /fell > 


(4.18) 


These spectra are sketched in figure [T] 

Since we have fixed the value of the spectral exponent a and since the spectra of (p 
in the phase-mixing-dominated region, where this exponent applies, propagate to higher 
m’s, we now also have determined Em{k\\^ fcj_) for fey > y/mk'^j_: see equation (14.81) . 

Physically, the validity of the argument that led to the last set of results (equa¬ 
tion (I4.15|) onwards) hinges on our ability to produce phase-space spectra that are con¬ 
sistent with substantial cancellation of the phase-mixing flux at fey < k’^ and thus with 
the majority of the free energy residing in the low Hermite moments. Note that it is 
actually not controversial that the phase mixing should be negligible for A:y <C k"^ be¬ 
cause this means the phase-mixing rate is low compared to the nonlinear advection rate, 
fcyi’th ^ k±u±, but, as we saw above, the energy is substantially dominated by the 
critical-balance curve fcy ~ where the two rates are comparable (recall our critique 
of the “fluid” theory in section r2.4.4D . In what follows, we shall build a case for the 
spectra that we have just derived—and so we will carefully avoid using the constant-flux 
argument (I4.15P and keep all the scaling exponents general. 


4.3. Spectra of higher moments in the advection-dominated region 
Let us now consider the higher Hermite moments, m ^ 1. The condition for the phase¬ 
mixing rate to be negligible compared to the nonlinear advection rate is the opposite of 
the condition (14.51) : 

k\\Vt\i^Vmk±u± fcy <C (4.19) 

(i.e., the part of the phase space on the side of the phase-mixing threshold opposite to the 
phase-mixing-dominated region). In section [4.11 the phase-mixing threshold was derived 
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by arguing that it represented the value of s up to which the “+” perturbations at some 
given k\\ and k± could propagate before being diverted to different wave numbers by the 
nonlinear coupling. More generally and more formally, we may simply argue that, if, say, 
the s dependence of / is a power law (which it must be; see appendix ICl) . the size of 
the phase-mixing term in equation (13.121) can be estimated as ~ (A:||Uth/s)/, which is 
negligible compared to the nonlinear term if the condition (14.191) is satisfied. 

The advection-dominated region 

fc||nth < k±u± ^ fc|| < k^j_, (4.20) 

which, as we argued in section 14.21 contains most of the energy content of the low- 
m Hermite moments (collectively represented by ip) and, therefore, of the advecting 
velocity uj^, is well within the domain of validity of the condition (14.191) . provided m ^ 1. 
Therefore, if we restrict our attention to the wave numbers (14.201) . we may neglect the 
phase-mixing term (second on the left-hand side) in equation (j3.12|l and thus deal with 
what is a purely “fluid” equation for /(s) at each s. As we argued at the end of section l5751 
we can then have solutions satisfying/*(fc||) = /(—/cy), for which the “-I-” and ” spectra 
are the same and the Hermite flux is zero (see equation (13.271) 1: 

E+m^,k^)^E-{kn,k^). (4.21) 

Physically, this is because, in the advection-dominated region, the nonlinear coupling 
between positive and negative fey mediated by the velocity field uj_ will be vigorous and 
fast—assuming, importantly, that interactions between uj_ and / are local in fey and so, 
in the right-hand side of equation (13.121) . the sum is dominated by wave-number 
triads with py ~ fey ~ fey — py. 

With the phase mixing neglected, the variance of / is (approximately) conserved at 
each s. The field /(s) is nonlinearly cascaded to smaller scales (larger k±) by the advecting 
velocity uj_, so the standard constant-flux argument gives us 

/^(s)fej_u_L ^ const(s) /^(s) oc fcj” => E:^{kj_) oc fej””^, (4.22) 

where E:^{k±) is the ID perpendicular spectrum of gm- Note that the spectrum has an 
m dependence, which cannot be determined via this argument. 

Since, in the advection-dominated region (j4.20p . the parallel-communication times are 
long compared to the nonlinear-decorrelation times, the perturbations can be expected to 
have a white-noise spectrum in fey. Therefore, we can write their 2D spectrum as follows 

^Ofe-d 

A™(fcy,fci)^A±(fey,fe^)^ fey<fel, (4.23) 

where we have allowed for an as yet unknown scaling with m. The fej_-scaling exponent 
d can be determined via the requirement that equation (14.231) be consistent with equa¬ 
tion (j4^2|) : assuming that, for any given fe_L, the region fey < contains most of the 
energylj 

feOfe“'^ 

dfey A„(fey,fe_L)/ dfey '' OC fe"^*^"”^ ^ d = l + 2r. (4.24) 

Since, as we have argued here, the Hermite flux is (approximately) zero in this region 
of wave-number space for all higher m’s, there can be very little net free-energy flow out 

t Technically speaking, we do not yet know this. We will justify this assumption a posteriori 
in section [4.4. II 
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of low m’s—and so, in retrospect, we were justified in assuming in section [4.21 that the 
energy of the low m’s was conserved and so a constant-flux argument Ij4.15|) could be used 
to deduce the scaling of ip. This allows us to adopt the scaling exponents (j4.16p (which 
we have thus far avoided using), in particular, r = 4/3, and so, using equation (j4.24p . 
the free-energy spectrum is (see figure [2] for illustration) 

E^{h,kx) ^ , fell < (4.25) 

Unsurprisingly, there is continuity between the scalings of and the scaling of Ei^. The 
scaling exponent a will be found in section H.4.11 

Before moving on to complete our scaling theory, we note that jumping to the result 
(|4.25|) already in this section was borne of pure impatience: we will discover in section HlTl 
that, in fact, it is possible to determine the scaling exponent d without relying on the 
as yet perhaps somewhat unconvincing claim that a constant-flux argument is legitimate 
for p despite the phase mixing being notionally not small along the critical-balance 
curve fc|| ^ 


4.4. Intermediate region and matching conditions 

We now have the form of the free-energy spectra in two regions, k\\ > y/m (phase¬ 
mixing dominated, very little free energy, equation (14.81) 1 and fey < (advection- 
dominated cascade, contains most of the free energy, equation (14.251) 1. It remains to 
determine the free-energy spectrum in the intermediate region between these two: 

U+(fe||,fc^)^ (4.26) 

(the following argument will only apply to the modes; both the reasons for this and 
the way to determine the spectrum of the ” modes will be explained in section EM- 
We have three new scaling exponents, but we also have the requirement to match 
equation (14.261) with equations (14.81) and (14.231) along the boundaries of the intermediate 
region. This gives us four relations 

a' = CT, d' -\- a'r = d, d' + a'r = ar — 3, a' + 2a' = 1 -|- a, (4.27) 


which we rearrange so: 


d-d' 


d = ar — 3, 


a = a = 


3 +r + d' 
'Yr ■ 


(4.28) 


Let us combine these with two equally uncontroversial (i.e., requiring no leaps of physi¬ 
cal intuition) matching and consistency relations from section EH] equations (14.111) and 
(I4.13p . which we can rewrite as 


8 -r 



The second of equations (14.281) then gives 


c = 5 — r. 


(4.29) 


d = 5-r = c. (4.30) 

Thus, the perpendicular scalings of Ei^ and E^ must be the same in the advection- 
dominated region k\\ < /c/). If we bring in equation (14.241) . i.e., the constant-flux argument 
(I4.22p . we get immediately 



a = 5. 


(4.31) 
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These are the same as the exponents (j4.16|) , except the need for the “fluid” constant-flux 
argument (|4.15|1 for ip has now been obviated by the combination of a more solid phase- 
space argument (|4.22p and a number of inevitable consistency relations. Equations (I4.28P 
will give us o', a' and a if we know d'. In order to determine the latter, we must now 
consider why physics in the intermediate phase-space region < fey < should be 

at all different from what happens in the advection-dominated region fey < considered 
in section 14.31 (and so why d' ^ d). 

4.4.1. Spectra in the intermediate region 

Except at the phase-mixing threshold fcy ~ i/m k'^^ (which is considered with more 
care in appendix [C|) , phase mixing in the intermediate region is dominated by nonlinear 
advection. However, interactions between uj^ and / cannot, unlike in the advection- 
dominated region discussed in section HTSl be local in fcy. Indeed, we know from section H?^ 
that there is very little energy left in uj_ at fcy ^ fc(J_. Assuming interactions to be local in 
fc_L, the wave-number sum in the right-hand side of equation (13.121) will be dominated by 
py < fc([_ (see appendix |B] for a careful analysis of the possible nonlocal interactions in the 
intermediate region). If fcy ^ fc(j_,thenpy -C fcy, /(fcy—py) ~ /(^y)j and so equation (I3.12|) 
now describes the advection of /(s,fcy) by an essentially two-dimensional velocity field 
(its parallel scale is much longer than that of /), with both s and fcy appearing as implicit 
parameters. This means that the variance of / will be conserved for each individual s 
and fcyU —and this in turn, by yet another constant-flux-in-fcj_ argument, implies 

/^(s, fcy)fcj_w_Lconst(s, fcy) ^ /^(s, fcy) oc fc/” ^ Em(fcy, fcj.) oc fc/’'^^ 

(4.32) 

This scaling is of the 2D spectrum, not of the ID perpendicular one, because fcy is a 
fixed parameter, rather than a variable over which there can be any nonlinear coupling. 
Comparing equations (14.321) and (14.261) . we read off d' and hence, with the aid of the first 
and third equations (14.281) . complete the determination of all scaling exponents: 

d' = r-|-l = ^, a' = l, a' = a = (4.33) 

O Z 

Thus, the free-energy spectrum in the intermediate region is 

^11 \/mfc^'^^, (4.34) 

modes will be discussed in section H31 

4.4.2. ID spectra 

In the run up to equation (I4.24L we assumed that integrating the 2D spectrum 
£’m(fcy,fc_L) with respect to fcy over the advection-dominated region fcy < k\^ captures 
most of the free energy contained in any fixed k±. Now in possession of equation (14.341) . 
we see that this is not entirely correct: in fact, using now both equations (14.251) and 

f This, incidentally, addresses a possible objection to the arguments in section HTSl that might 
have been troubling a perceptive reader: were we really justified in assuming that the conserved 
variance of /(s) could all be accounted for within the region k\\ < fc/, leading to the constant-flux 
argument (14.221) ? The answer is that the free energy outside that region is either subdominant 
(section |4.1l) or conserved separately (section 14.4.11) . 




fc-lfc-'^/^ 

fin 
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5/2 


sketched in figure [51 The spectrum of the 
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(I4.34L we find the ID perpendicular spectrum to be 
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E^'^{kj_) = j dkiiE+{kii,kj_) 


i.4/3 


k 


m 3/2 7.4/3 


u4/3 


d/cii 




m 


5/2 


/-7/3 


m5i2 


(1 + . 

(4.35) 


So there is logarithmically more free energy in the intermediate region, but this does not 
affect the k± scaling, which is what we were after in equation (I4.24|) . so the derivation in 
section 14.31 survives. 

For completeness, let us also calculate the ID parallel spectrum. Integration over k± is 
dominated by the wave numbers around the phase-mixing threshold k± ~ (fc|| fy/rn)^^^, so 


E^Jrt(.h) = J dkj_E+{ki\,kj_) 


j 


dk± 


kik-^ 


1 , 3/4 


/(fcii/^A/;)3/4 


dfc_L 


^11 ^-L 

777,3/2 


(4.36) 




Note that this m ^ scaling appears to be in decent agreement with the Hermite-space 
spectra reported bv iHatch et al\ ( 20131 1201^ . 


A perceptive reader might be feeling a growing resentment over our use of the spectrum 
(|4.34l) at wave numbers around the phase-mixing threshold fey ~ y/mk'^^, even though, 
technically speaking, we have only justified equation (j4.34|) in the region k'^^ < fey <C 

y/mk'^^, seeing that at fcy ^ y/mk'^^, phase mixing cannot be neglected compared to 
the nonlinear advection. In appendix [Cl we show that it nevertheless makes sense simply 
to match the spectrum (|4.26|) to the phase-mixing-dominated spectrum (14.8|) along the 
phase-mixing threshold. 


4.5. Anti-phase-mixing spectra 

The arguments about the intermediate-region spectra presented in section 14.41 only ap¬ 
ply to the spectrum of the “-I-” modes. Since the advection velocity is effectively two- 
dimensional in the intermediate region (see section [4.4.IL there is no coupling between 
different parallel wave numbers and so no echo effect. Thus, if, as we argued in sec¬ 
tion 231 Ri £1+ oc in the advection-dominated region, fc_L > '^'^e to vigorous 

coupling between parallel wave numbers, and in the intermediate and phase- 

mixing-dominated region, fcj_ ^ fcy , due to absence of any such coupling, we must 
expect that the energy-containing wave numbers for the ” modes are ones along the 
critical-balance curve fc_L k^\'~ ■ 

The fc_L —>■ 0 asymptotic behaviour of the ” spectrum must be the same as for any 
other field, oc because the reasons for it are purely kinematic fappendix El). 
Thus, we posit 

fc3 

A-(fcy,fc^)~^i^, fcy>fcl (4.37) 

and impose matching conditions between this spectrum and equation (14.231) at fcy ^ 


a" = ^ = 5, 


a a 


(4.38) 


r 
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Figure 2. Spectra of higher Hermite moments, i?^(fc|| , k±)-. (a) in the (fc||, m) plane at constant 
k± {E^ and are shown in the right and left panels, respectively; in the left panel, fey increases 
leftwards), (b) in the (fc_L,fc||) plane at constant m {E^ and E^ are shown in the upper and 
lower panels, respectively; in the lower panel, fey increases downwards), (c) vs. fcy at constant k± 

and m, (d) vs. k± at constant fcy and m, (e) vs. m at constant fey and k± (such that fcy > k'^^ , 
otherwise the spectra are E^ ^ at all m). All plots are logarithmic. The spectrum 

i5,p(fcy,fcx) (see figure[TJ is given in (c) and (d) as a dashed line, for reference. Free-energy flow 
through phase space as represented in (a) and (b) is described in section [53] 








































Phase mixing vs. nonlinear advection in plasma turbulence 


29 


This completes the determination of the phase-space spectra of the anti-phase-mixing 
component of the free energy: 


E^{kii,k±) 


1 




m 


5/2 


k?.k 


- 11/3 


if fc|| > k 
if fell 5- k 


4/3 
_L ’ 
4/3 


(4.39) 


Figure [5] shows these and illustrates their relationship to other spectra derived above. 
The ID spectra that follow from equation (14.391) are 


<ik±E-{kii,kj_) 


k 


- 7/3 




m 


5/2 ■ 


(4.40) 


(4.41) 


Note that these are both subdominant, in m, to the “-(-’’-mode spectra ()4.35p and (14.3611 . 

4.6. Effect of collisions 

4.6.1. Collisional cutoff for phase-mixing modes 

In the phase-mixing-dominated regime 1section [4.1|) . the collisional cutoff is set, in the 
same way as in the linear theory (section 12.Sp , by the competition between the phase¬ 
mixing rate ^ k^^Vth/V^ the collision rate ~ ixm. The perturbations are collisionally 
damped if 


> 


k\\Vth 


> k±u± 


> 


2/3 


k^ 

fell > 


(4.42) 


giving a cutoff in Hermite space (cf. equation (I2.59l) ill1 In both the intermediate (sec¬ 
tion |T4|) and advection-dominated ('section IT^ regimes, the relevant comparison is be¬ 
tween the collision rate and the nonlinear-advection rate: 

l4/3 


vm > kxu±_ P 


^11 




k 

> ^ 


< 


(4.43) 


These cutoffs are sketched in figure O 

At a fixed m, the above relations imply that there is an infrared collisional cutoff in 
the (/cy, k±) space: perturbations are damped if 


fc|| < 12771^/^, k± < (yrrifl^. 


(4.44) 


These cutoffs will not, of course, be relevant in comparison with the outer scales (/cyo and 
k±o', see section [2.4.II) except at high enough m or if the collision frequency approaches 
the characteristic phase-mixing and nonlinear-advection rates at the outer scale. In the 
latter case, one expects some amount of free energy to drain via collisional dissipation 
around the oute r scale, petering ou t at larger (fe y, fc i ), where the collisionless physics 
asserts itself (cf. [Hatch 201Ial f6. 2013 . 2014 ). 


4.6.2. Collisional cutoff for anti-phase-mixing modes 

The collisional cutoff (14.431) on the free-energy spectrum of the “-I-” modes in the 
advection-dominated regime must extend to the ” (anti-phase-mixing) modes because 

f In the last formula in equation (14.4211 . we implicitly nondimensionalised the collision fre¬ 
quency: i'Z/||/uth = h||/Amfp —> n, so the Hermite cutoff is me ^ (fey Amfp)^^^, where Amfp = vth/n 
is the mean free path. In other words, in rescaled units, one can replace v 44- 1/An,fp wherever 
this makes things more transparent. 
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Figure 3. Partition of phase space, viz., (fc||,m) plane at fixed k±, showing the 
collision-dominated region. Cf. hgure[2|a). All axes are logarithmic (in the left panel, fey 

increases leftwards). 


nonlinear coupling is the only source of the latter (see section H751) . But the anti-phase- 
mixing modes propagate from high to low m and so a zero “boundary condition” at 
high m will be imprinted onto a region of phase space at lower m’s. To wit, arguing 
analogously to section HIT] and considering now fc|| < 0 in equation (13.121) . we note that 
anti-phase-mixing modes propagate along the characteristics 


s 


l^llkth 

V2 


t + So, 


(4.45) 


where so is a constant. Whatever anti-phase-mixing spectrum exists at s = so, it will 
be replicated over all s satisfying equation (14.451) for times shorter than the nonlinear 
time, t ^ T^\ ik±u±)~^ ^ Assuming that is cut off for s > Sc = k^^ j^/v 

(equation (14.431) ) and letting so = Sc, we conclude that A” must also be cut off for 


_ \k\\\vth 

k±u± 






lfclll 


rUc 



(4.46) 


This implies, in particular, that there is no anti-phase-mixing energy at any m’s for wave 
numbers satisfying 




(4.47) 


The collision-dominated region of the phase space is sketched in figure El 
Whereas the collisional cutoff is safely removed to infinite m in the limit u —> 0, in 
systems with only moderately low collision frequency, one should expect to see a finite 
reduction in the anti-phase-mixing flux at higher m’s, as per equation (I4.46p . 


4.7. Total free energy and dissipation 

In the linear problem, where all energy injected into the system had to be removed by 
Landau damping (meaning phase mixing followed by collisional dissipation at high m), 
the free energy stored in phase space in a steady state had to diverge with vanishing 
collisionality (see equation (I2.61|) l in order for the dissipation to remain finite (equa¬ 
tion (12.601) 1. In the nonlinear situation with which we are now preoccupied, the Hermite 
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spectra are steep power laws and so the free energy will be finite and collisional dissipa¬ 
tion vanish, with all of the injected energy having to be removed via dissipation at small 
spatial scales (sub-Larmor and so outside the regime of validity of this theory). 

To demonstrate this a little more quantitatively, let us repeat the calculation of the ID 
parallel spectrum (equation (14.361) '). but now, in integrating the 2D spectrum over k±, we 
assume that there is no free energy at perpendicular wave numbers below the collisional 
cutoff kj_ ~ (vrnf'^'^ (equation (14.441) '). The following three cases correspond to the 
collisional cutoff falling into the phase-mixing-dominated, intermediate and advection- 
dominated regions, respectively: 




Note that the last two cases are only relevant at very high to (because fey ^ feyo ~ 1, the 
outer scale in our units). Now integrating these spectra over to, we find that the total 
free energy in a given fey is completely dominated by low to’s: 


if TO < [ — I ~ TOf, 


■-1 


vrrP 

1,0 

^11 


if 

V j V 


if TO>^. 

V 


(4.48) 


/ oo 

^ dTOi;i+(fcy) - 

The total collisional-dissipation rate vanishes with v. 


(4.49) 


/'™° dm /fell \ 

dmmEj^(k^^) ^ k^'^v / -~ In f — j —^0 as —>■ 0. 

(4.50) 

Equations (14.491) and (I4.50p are the nonlinear versions of equations (12.611) and (12.601) . 
respectively—we see that, unlike in the linear problem, the free energy remains finite 
and collisional dissipation vanishes as z/ —>■ -|-0. 


5. Conclusion 

5.1. Summary of free-energy spectra and of the method of deriving them 

Considering the full phase space (fcy, fc_L, to), we posited a set of power-law relationships 
for the free-energy spectra and then determined the scaling exponents from a combina¬ 
tion of matching conditions between different regions of the phase space and physical 
arguments about the free-energy flows, constrained by conservation laws. The spectra 
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are, for the phase-mixing modes (i.e., perturbations that propagate from low to high m), 




^11 ^-L 

^ 5/2 

ku^k^, 


if fell 5-, fc 


,4/3 


(advection dominated), 

^ ^11 ^ \/m k'^^ 
(intermediate), 


(5.1) 


if ki\ > 


m 


^4/3 


(phase-mixing dominated), 


for the anti-phase-mixing modes (propagating from high to low m), 


E- 


, fcj.) 


J.0j,-ll/3 

™ 5/2 


^-5iL3 


if fcii < fe 


,4/3 


(advection dominated). 


(5.2) 


m 


5/2 


if fell k 


4/3 


(no echo), 


and for the “fluid” (low-m) moments, 


if ku < k^/^ 


, k±) 


(advection dominated). 




if k\\ > k 


4/3 


(5.3) 


(phase-mixing dominated). 


A graphical summary of these spectra is presented in figure [2] 

As is manifest in the above formulae, the phase space is partitioned into several regions, 
where different physics controls the distribution of the free energy. 

• In the phase-mixing-dominated region lsection l4.lL the phase-mixing rate is greater 

than the rate of nonlinear advection, k\\Vt\i /^ k±^u_\_, and so whatever distribution 
of free energy exists at these wave numbers at low m’s will simply be propagated to larger 
m’s—this is the part of the wave-number space where modes are “Landau-damped” in the 
usual linear sense. The perpendicular spectrum in this region (oc see equation (15.11) 1 
is fixed on purely kinematic grounds (appendix]^, the m scaling (oc is the same 

as in the linear problem, c orresponding to constant Hermite flux (IZocco fc Schekochihin 
2011 . Kanekar et a/.|[2015l : see section [231) . whereas the scaling exponent of the parallel 
spectrum (oc a = 5) is hxed by matching with the nonlinear dynamics Isection 14.41 
see equation (14.311) 1. 

• The continual flow of free energy into high m’s as described above sets the matching 
condition at the phase-mixing threshold, where the nonlinear advection rate becomes 
comparable to the phase-mixing rate, k±u± ~ k^^Vth/V^- The role of the nonlinear 
advection is to divert the free energy from flowing straight to higher m’s to flowing to 
higher k±’s. The competition between these two processes sets the prevailing dependence 
of the free energy on m, giving rise to the m“®/^ scaling of the 2D spectra (equation (15.IL 
derived in section 14.4.11 and appendix [C]) and the m“^ o verall scaling (equation (I4.36L 
in reasonable agreement with recent numerical studies by Hatch et '^ 120131 [2nil . The 
situation at the phase-mixing threshold is so crucial because the free-energy spectra rise 
as fc_L increases and fcii decreases from the phase-mixing-dominated region towards the 
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phase-mixing threshold and then fall beyond it, at higher kj_ and lower fey, so it is along 
the phase-mixing threshold that the energy-containing scales in phase space lie. 

• The intermediate region comprises the wave numbers at which the nonlinear-advection 

rate is already dominant compared to the phase-mixing rate of the high-m moments of 
the distribution function, but not of the low-m moments and, in particular, of the zeroth 
moment, (^, which is what sets the E x B flow velocity Uj^ that is doing the advection: 
k\\Vi\i/y/m fc_LU_L <g; The energy-containing wave numbers for the flow lie along 

the critical-balance curve k±^ux ^ k\\Vt\i Isection 14.21) — and so the nonlinear interactions 
in the intermediate region are nonlocal in /cy, with short-parallel-scale perturbations of 
the distribution function advected by a longer-scale flow, i.e., an effectively 2D veloc¬ 
ity field (see appendix |B|) . A constant-flux argument for the free-energy cascade in k±^ 

_ '7 /O _-1 

then fixes the fcj_ scaling of the 2D free-energy spectrum in this region, while its fey 
scaling follows from matching to the spectra at the phase-mixing threshold and at the 
critical-balance curve (the second scaling in equation (15.11) : see section lT4l) . 

• Beyond the critical-balance curve, ^ fcycthj the nonlinear advection is com¬ 

pletely dominant over phase mixing, giving rise to the advection-dominated region. The 
advecting flow is now 3D and another constant-flux argument gives the scaling of 

the free-energy spectrum, whereas its k?, scaling is a white-noise spectrum deduced via 
a simple causality argument implying that perturbations with a certain perpendicular 
scale are decorrelated at parallel distances long enough that information cannot traverse 
them at the speed ^ Uth over one cascade time corresponding to that perpendicular scale 
(the first scaling in equation (15.11) : see section IT31) . 

The above arguments have all focused on the free energy contained in the perturba¬ 
tions that propagate from low to high m, i.e., ones prone to phase mixing (whether it 
is fast or slow compared to nonlinear advection). In a nonlinear system, an advecting 
flow that has a parallel spatial dependence, i.e., fcy ^ 0, can couple these perturbations 
to others that have parallel wave numbers of opposite sign and so will propagate from 
high to low m, a phenomenon known as plasma echo Isection 13.21) . Separating all per¬ 
turbations into these “-I-” and ” components (section 13.11) allows us to express the 
free energy as the sum of their spectra and its flux in Hermite space as proportional to 
the difference between these spectra (section [3^. In the advection-dominated region, 
vigorous nonlinear coupling implies that the “-I-” and ” spectra are the same and so, 
statistically, there is no free-energy flux between different m’s—i.e., the phase-mixing and 
the anti-phase-mixing energy fluxes cancel each other (see section 14.31 and the first scal¬ 
ing in equation ()5.2I) 1. In contrast, there is no echo effect and, therefore, no significant 

” energy either in the intermediate region (because the flow velocity there is effec¬ 
tively 2D and so cannot couple different fcy’s) or in the phase-mixing-dominated region 
(because anti-phase-mixing modes do not propagate to higher m’s). The ” spectrum 
outside the advection-dominated region (the second scaling in equation (15.21) . derived in 
section 14751) is, therefore, determined by the kinematic constraint giving the k\ scaling at 
long wavelengths and by the matching conditions along the boundary of that region—the 
critical-balance curve. 


Finally, the spectra (15.3|) of the low-m, “fluid” moments are basically a continuation 
of the high-m spectra dsn and (|5.2|) down to low m’s. Physically, since the Hermite 
flux between different m’s is on average shut down in the advection-dominated region, 
these scalings can be determined by assuming constant flux of the “fluid ” part of the 
free-e nergy, i.e., effectively, by pretending that the turbulence is fluid-like (iBarnes et al 


2011 


see sections [2. 4. 2l and l4.2|) . Such a shortcut has always been tempting le.g.. lWeiland 


19921) . but was not a priori justihed for a kinetic system (section [27474]) . 
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5.2. Free-energy flows 

Although this was implicit in our discussion of the partition of phase space ('section f5.ll) . 
it is worth spelling out what path the free energy takes through it. Let us start from 
some (m, fc||,fej_) in the phase-mixing-dominated region (lower right corner in the right 
panel of figure [SJa)). At first, the free energy will move (phase-mix) from there to higher 
m (vertically towards the blue line in figure UKa)) until it reaches the phase-mixing 
threshold m ~ (equation (14.5(1 with r = 4/3; the blue line in figures [2Ka,b)). 

There it enters the intermediate region, where it is advected by an effectively 2D velocity 
field (see section 14.4.11 and appendix |B|) to higher k± while staying at fixed fc|| (in the 
top panel of figure [2Kb), horizontally from the blue towards the red line) until it reaches 
the critical-balance threshold /cj_ ~ (equation (I4.20p : the red line in hgureUKb)). At 
that point it enters the advection-dominated region, where the advection is 3D and the 
energy flows along the critical-balance curve (diagonally upwards along the red line in 
the top panel of figure |2Kb); see section |43)) . Since a 3D velocity is effective at coupling 
positive and negative kfs, this flow of energy involves both “-I-” and ” modes (the 
latter shown in the bottom panel of figure |2Kb), where the critically balanced energy 
flow is also along the red line). There is not much flow of the ” energy beyond the 
critical-balance threshold (to the left of the red line in the bottom panel of figure HKb), 
or, equivalently, to the lower left of the red line in the left panel of figure |2Ka)) because it 
nonlinearly couples back to “-I-” modes faster than it can anti-phase-mix to lower m’sl3 


5.3. Implications and outlook 

The free-energy distribution in phase space summarised above has several important 
properties and implications. 

The free-energy flux out of the “fluid” moments is heavily suppressed in the wave- 
number region bounded by the critical-balance curve, /cy < k'^^, which is also the region 
that contains most of the free energy flowing through the inertial range. Thus, at the 
energetically relevant wave numbers of the inertial range, Landau damping is effectively 
absent. The resulting Hermite spectra have steep power laws (oc for the total energy; 
see equation (14.361) 1 and so the total free energy contained in the phase space is hnite, 
dominated by low m’s (the energy in the “fluid” moments) and does not diverge at 
vanishing collisionality (equation (14.491) 1 —in sharp contrast to its behaviour in the linear 
problem (see equation (I2.6ip l. Furthermore, the total collisional dissipation vanishes in 
the nonlinear problem (equation (14.501) 1. again in contrast to the linear case, where the 
dissipation rate is hn ite and absorbs all of the energy that is injected into the system 
( Kanekar et al ][2H5t see equation (12.601) 1. This means that most of the dissipation occurs 
at small spatial scales (i.e., beyond the Larmor scale, a region that we have left outside our 
detai led focus). This is indeed what was recently found numerically bv [Hatch et al. ( 2013 . 
20141) : decreasing share of the collisional dissipation with decreasing collisionality. Note 


that hnite collisionality imposes a cutoff on the free-energy spectra at high enough m’s, 
or, equivalently, at low enough fey and k± (equation (14.441) 1: when the collision frequency 
approaches the rates of phase-mixing and nonlinear-advect ion rates, a certain amount 
of collisional dissip ation will occur at low wave numbers (cf. IWatanabe fc Suefamall20(i5 
Hatch et all 2nilfl 3). 


It is inevitable that one must ask about the implications our results might have 


t It is possible for mode-coupling in fey to combine with anti-phase-mixing to push some 
” energy towards larger fey, but that process is diffusive in fcy and will be slower than direct 
nonlinear coupling back into modes. It becomes important w hen the advecting flow i s 
scale-separated from the distribution function that is advected by it (ISchekochihin et al.ll2ni(i) . 
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for the Landau-fluid closures as a viable modelling techniqu e—a subject that has long 
been discussed and refined in the context of fusion plasmas dHammett fc Perkins 1^90 1 


Hamrnett et al.lll992 . 19931; Weilandlll992 : Mattodll99^ Hedrick fc Leboeullll992t Dorland fc Hammett 

1993: Beer fc Hammetd 1996t Snvder et^ 19971: Snvder fc Hammett 200 loll 6: Ramos 

2005 ) as well as. more recently, space and astro physical ones ( Passot fc Sulemll2004 . 20061 


2007 : Goswami et al. 20051: Passot et aZ.ll2012l) . While the basic idea of the Landau-fluid 
approach is to include into fluid equations damping terms |/c|||uth) fit to capture cor¬ 
rectly the linear Landau damping, it has long been known in this field that quantitatively 
these models work better when more Hermi te moments are retained and this inclusion 
happens at the level of the highest of them I Smith 1997t) . Considering that the free en¬ 
ergy scales steeply with m, as shown above, it stands to reason that, at low collisionality, 
Landau-fluid closures that retain a certain finite (independent of the collision frequency) 
number of moments may be sufficient for a full characterisation of kinetic turbulence—in 
that already just this finite number of moments will be enough to capture most of the 
echo flux from phase space back to “fluid” moments. The Landau closure terms affecting 
the highest of the retained moments will then serve to regularise the problem in the ener¬ 
getically subdominant part of the wave-number space—the phase-mixing region—where 
the free energy has a shallow scaling ~ ('section l4.ll) . 

There is clearly space for further development of this line of reasoning, leading to 
more quantitative prescriptions for capturing the echo effect within the Landau-fluid 
framework. If one thinks of these closures in the same mo delling spirit as on e does 


about Large-Eddy-Simulation techniques in fluid dynamics (iSmagorinskv 


more recently, in gyrokinetics (iMorel et al.ll201lL l2012t iBanon Nayarro et al. 


—and, 


20141) .— the 


j.j.j-5/2 spectrum we haye deriyed can serve the useful role of providing the signature of 
a well-developed nonlinear phase-space “cascade,” which, once formed, can be promptly 
and safely cut off by model dissipation terms. 

Ranging somewhat further afield, we note that spacecraft measurements of compressive 


turbulence (ICelnikier et all 19831 1987: Marsch & Tul 1990l: 

Bershadskii & SreenivasanI 

I 2 OO 4 I: iHnat et al. 2005: Kelloe-ff & HorburvI 2005: Chen et al. 

2011. 2014) show healthy 

Kolmogorov-like power-law spectra—in what is generally a 
Landau damping of such fluctuations (iBarnesI 196fil) ought 

P 1 plasma, where the 

to be of the same order 
imilarlv robust power-law 

as their nonlinear cascade rates (Schekochihin et aZ.ll2009r). 8 

spectra at sub-ion-Larmor scales have also been measured (f 

iahraoui et al. 20091 2010t 

2013: Alexandrova et all 20091 20121: Chen et al. 2010t 2013) 

and found in kinetic sim- 
h Landau damping of ki- 

ulations (Howes et al. 2011: IChantr et al. 2011). even thoug 

netic Alfven waves (Howes et_a^l2006l Carv & Borovskv 200f 

1) should be quantitatively 

noticeable at these scales (Howes et al. 2008: Podesta et all 

2 OIOII. Whereas attempts 

have been made to argue that in some of these situations the linear damping might be 

weak (Lithwick & Goldreich 200ll iHowes et al. 20081: ISchekochihin et al. l2009l). it is a 


tempting—and more interesting—thought that the general mechanism for (statistical) 
suppression of phase mixing in a turbulent system proposed here is responsible for making 
collisionless plasma turbulence in the solar wind behave in a seemingly more “fluid-like” 


t These can be shown to be drift-kinetic fields passively advected by the turbulent velocity 
field ux associated with Alfvenic perturbations. They satisfy equations that are quite similar to 
equation (I2.8|l . although with an additional complication that particles stream along magnetic 
field that is also pertu rbed by the Alfvenic turbulenc e and so linear an d nonlinear mixing are 
somewhat intertwined (|Schekochihin et al.ll2007l . 120091 : feunz et al.ll2015l 'l. 
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fashion than theoreticians might have thought it had a right to doi A numerical and 
theoretical investigation of this possibility is a subject of our current efforts. 

To conclude, the considerations presented above appear to point to a number of promis¬ 
ing directions for numerical experiment and further thought. We hope to explore some 
of those in the not so distant future^ 
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DoE grants DE-FG02-93ER54197 and DE-FG02-08ER54964. All authors are grateful to 
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Appendix A. Long-wavelength scaling of spectra 

Here we review the standard argument that the spectrum of a 2D-isotropic homoge¬ 
neous field (^(fc||,rj_) has the low-wavelength asymptotic form 

^.l) = 27r/c_L(|(/5(/c||, k_L)p) oc as fcj, —>■ 0. (Al) 

Setting the perpendicular-Pourier-transform conventions to be 

V?(A:||,k_L) = y 

where L± is the box size, we find 

[ drj_rj_Jo{kj_rj_)C{k\\,rj_), (A4) 
^_L Jq 

where rj_ = rj_i — rj_ 2 , C(fc||,r_L) is the two-point correlation function of ip (which only 
depends on rj_ = |rj_| because the field is statistically homogeneous and isotropic), and 
Jo is the Bessel function of order zero. 

j: An immediate physically interesting conclusion from such an outcome, apart from power-law 
compressive spectra being theoretically legitimised, would be that one should not expect any 
ion heating associated with the inertial-range turbulence (see equation (14.501) 1. with the thermal 
fate of all turbulent energy determined at the ion Larmor scale, where t he 4D drift-kinetic 
phase-space cascade morphs into a more complicated 5D gyrokinetic one ([Schekochihin et al\ 
120091: iHowes et al.ll201ll: iTold et al.llMll . 

t As this p aper is going int o press, the firs t de dicated numer i cal te sts of our theory have been 
undertaken bv iKanekarl (120151 1. [Parkerl (120161 1 and iParker et al.1 (120161 1. so far broadly supporting 
our conclusions. 
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If the correlation function C(fe||,rj_) decays sufficiently quickly with rj_, it will restrict 
the integral in equation (lA 4|) to values of r_L that are smaller than or comparable to 
the perpendicular correlation length r±c of the field (/3(fc||, rj^). Note that rj_c will be a 
function of fey, so it is not necessarily the outer scale—in section 021 we argue that it is 
the critical-balance scale, rj_c ^ fey ■ If we now consider kj_rj_c ^ Ij we may 

expand the Bessel function Jo{k±r±) = 1 — /A + ... in equation (IA4I) . which then 
gives us 


E^{k\\,kj_) = 27r/c_L(|(p(/cy,k_L)|^) = k_i 


+ k 


3 

_L 



dr_Lr_LC'(fcy,r_L) 


dr±r]_C{ki\,r±) 


(AS) 


The first term is proportional to J d^rj_((/5(fcy,rj^)(/j*(fcy,0)) and so it vanishes if we 
assume that f d^rj_‘/?(fcy, rj_) = ip(kji,k± = 0) = 0, i.e., that there are no purely ID 
parallel modesl3 Hence we obtain the desired result (lA II) . with the proviso that C{k \\, ru) 
decays faster than 1 /r^ as rj_ —>■ oo and so the integral prefactor of k\ in equation (lA 5|) 
converges. 


Appendix B. Nonlocal interactions in the intermediate region 

Consider the nonlinear coupling expressed by the right-hand side of equation (I3.12p . 
which we now rewrite as a wave-number convolution in both parallel and perpendicular 
directions: 

()f\ 

— j (fcy,k_L) = -ik_L • ^ Uj.(py,p_L)/(A:y -py,k_L - pu). (B 1) 

/ nl PlhP-L 

In the intermediate wave-number region between the phase-mixing threshold and the 
critical balance, 

y/mk±u± > fcyUth ^ k±u± ^ y/rnk'^ > fcy > k'^, (B 2) 

the coupling in equation (EU must be predominantly between disparate wave numbers 
(i.e., the coupling is nonlocal) because the energy-containing wave numbers for uj^ are 
py < p(l, which lie outside the region (IB 21) . There are two basic possibilities: coupling 
that is local in k± but nonlocal in fcy and coupling that is local in fcy but nonlocal in k±. 
In analysing the rates of such interactions, we will consider / to be at the phase-mixing 
threshold, fcy ~ and uj^ in critical balance, py ~ p(j_. 

Suppose the perpendicular coupling is local, p_L ~ |kj_ — pj_| ~ k±. Then 

py ~ /cl ~ i < fcy, (B3) 

" y/m 

so the distribution function /(fcy ~ P||) ~ /(^||) is advected by an effectively two- 

t In the theory of a passive scalar, the quantity d^rxC(fc||, r^) is known as the lCorrsinI 
(119511 1 invariant—the decay laws for a passive scalar can depend on whether this invariant is zero 
or fin ite because that sets the long-wavelength as ymptotic behaviour of the scalar’s spectrum 
fe.g.. lEvink fc Xlnll200(]l : ISchekochihin et al\ 200-3 1. The fact that this asymptotic behaviour is 
~ k\ in our theory, will have implications for the determination of the Hermite spectrum; see 
appendix [C] 
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dimensional velocity field: back in real space, equation (|B 1() becomes 

«-uj_(z = 0,r_L) • V_L/(A:||,r_L). (B4) 

nl 

The rate of nonlinear advection of /(fc||) is, as usual, 

k±u± ^ (B 5) 

Note that as there is no coupling in fey, there can be no echo. 

Now suppose instead that it is the parallel coupling that is local, p\\ ^ |fc|| — P||| fc||. 

Then 

Pj_ ^ ^ » fc_L, (B6) 

so the distribution function / is advected by a much-smaller-scale (in the perpendicular 
direction) velocity field. The net effect of such an advection will be turbulent diffusion 
of / with the effective mixing length ~ 1/p_l and the effective diffusion coefficient 

Aurb ~ pT" - fcr (B 7) 

P± 

The rate of nonlinear advection associated with this process is then 

Aurb^i ^ ^(2_t)/2r' 

provided r < 2 (which it is, considering it will turn out to be r = 4/3). This is much 
smaller than the local-in-Aj nonlocal-in-fc|| advection rate (IH3- Thus, the latter type of 
interactions will be the dominant ones—the claim we make in section 14.4.11 which this 
appendix is meant to back up. 

Note that other kinds of interaction—of various degree of non-locality in both /cy 
and k± —cannot prove faster because non-locality in k± will always slow down coupling 
(diffusion is slower than advection) while more or less non-locality in fcy simply makes 
the velocity uj^ more or less two-dimensional compared to /, without changing the rate 
of advection. 



Appendix C. Spectra near phase-mixing threshold and the 
free-energy decay in Hermite space 

In a statistical steady state, the free-energy spectrum is independent of time and so 
described by equation (13.261) : 

(r(fcy)u^(py) . V^/>y -py)) - 2z.sA. (C 1) 

Let us consider the wave numbers around the phase-mixing threshold, for which fcyUth 
fcj_u_L, so the phase-mixing term is comparable to the nonlinear term. Ignoring colli¬ 
sions, assuming locality in k± (see appendix E]) and expanding in py ~ k±u±/vth “C 
y/rnk±u±/vth ~ fcy, we have (cf. equation (IB dp i 

(l/>ll)l') ^ = 0,r^). Vi|/>y)|2) . (C2) 

Formally, this looks like an equation for the spectrum of a passive 2D field /(fcy,rj_), 
parametrised by fcy, advected by a 2D velocity field vn_{z = 0,rj_) and decaying with s, 









Phase mixing vs. nonlinear advection in plasma turbulence 


39 


which plays the role of time. While devising a specific quantitative closure for the triple 
correlator in the right-hand side of equation (jC 2() is outside the scope of this paper, it 
is plausible that the solutions around the phase-mixing threshold sk±u± will 

satisfy, roughly, 


els 


k±u± j2 ^ 
k\\Vth s 




(C3) 


Thus, the decay must be a power law, as we indeed assumed in equation (14.261) . This 
decay law is set at the “outer scale”, which is the phase-mixing threshold: kj_ ^ (fc||/s)^/” 
(fc|| is fixed). Below this scale, i.e., at k± ^ (fcii/s)^/”, the phase-mixing term is small 
and the /(s,/c||) is simply cascaded subject to the constant-flux argument proposed in 
section [4.4.11 (i.e., the right-hand side of equation (1C 21) must vanish to lowest order in 
1/s). This gives a spectrum of the form (j4.26l) . inheriting its decay law ^ 1/m'^ from 
the “outer scale”. The decay law of the spectrum in the advection-dominated region 
k± ^ ^ equation (I4.25L is then the same, a = a', via matching at the critical-balance 

curve fc_L ^ (see equation (14.281) 1. 

What is the relationship between /r and a' and how is this scaling exponent determined? 
Equation (1C 3p effectively sets the ID parallel-wavenumber spectrum, i.e., as explained 
above, the free-energy content of all wavenumbers A:j_ > (/cn/s)^/”: using equation (|4.26|) . 
we get 






/(fc||/\/rn)i/'- 


dfc_L£'+(fc||,fc_L) 




-.g' — {d' —1) /2r ’ 


(C4) 


so fj, = 2a' — 1 — {d' — l)/r. This decay exponent, or, equivalently, a', is deduced (along 
with a') by matching the decay law del with the decay law of the total variance of 
/^(/cy) contained at long wavelengths k± < (/cy/s)^/”: using the asymptotic form 
we get 






dkjLE+{k\\,kjL) 


j^l/2+2/r^ 


(C5) 


SO fj, = 4/r. Matching equations (IC4p and (1C 51) . we get two relations constraining a' and 
a , which, combined with matching conditions at the critical-balance curve k± ^ fc|| , 
are the same as equations (14.281) . Note that using the set of exponents (|4.31l) and (14.331) . 
we happily recover the ID parallel spectrum (I4.36P from either of equations (1C 41) and 
(|C 51) . We also Hnd that ^ = 3. 

Note that deducing the decay law of a turbulent field by fixing the long-wavelength 
asymptotic behaviour of its spectru m (Em. oc in ou r case) is a standard trick of 
the trade in turbulence theory (e.g., 


Evink fc Xin 20001: Schekochihin et dll 


1941 a : CorrsinlflQ^ ISaffmanIIl967l : 


2004 Davidson 


201(l[l2013ll. 
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